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CHAPTER 1 
THEORY 

1.1 Pcirticle Physics 

Particle physics is the branch of physics dedicated to the study of matter 
and energy at the most fundamental level. This means that the job of a particle 
physicist is to identify the smallest constituents of matter and describe how they 
interact with each other. 

Humankind has been interested in this subject since ancient times. Two of the 
first particle physicists in recorded history are the Greek thinkers Empedocles and 
Democritus from the fifth century BC. Empedocles stated that our complex world 
was made from combining four fundamental elements (earth, air, fire and water) in 
different proportions. Democritus on the other hand, believed that the apparently 
continuous objects in the natural world were not really continuous, but made from 
voids and indivisible particles called atoms. 

Over the last half century Particle physics has advanced tremendously and we 
now have a beautiful, but incomplete, theory firmly grounded on experiment that 
describes the fundamental constituents of matter and how they interact with each 
other. This theory is called "The Standard Model." 
1.1.1 The Standard Model 

According to the Standard Model, the fundamental building blocks of matter 
are point like particles which interact with each other in as many as three different 
ways. Each type of interaction, or force, is itself carried by point like particles 
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called force carriers. The particles which are force carriers are bosons ^ and are 
collectively called gauge bosons because they are needed for the theory to be gauge 
invariant. The non-force carrying particles are fermions. ^ 

The three interactions described by the Standard Model are called the 
electromagnetic force, the weak force, and the strong force. The force of gravity is 
not included in the Standard Model, and this is one of the reasons the Standard 
Model is not yet complete. The electromagnetic force carriers are photons, the 
weak force carriers are the W~^, and W~ particles, and finally, the strong force 
is carried by eight kinds of gluons. Table II II summarizes the situation. 

Table 1-1: Gauge bosons and the force they carry. 



Symbol 


Name 


Force Carried 


7 


Photon 


Electromagnetic 






Weak 




w+ 


Weak 


w- 


w- 


Weak 


9 


gluon 


Strong 



Particles which interact through a particular force are said to couple to it 
and to carry an associated charge. The nomenclature is as follows, particles that 
interact through the electromagnetic force have an electromagnetic charge called 
electric charge, those which interact through the weak force have a weak charge 
called weak isospin, and those which interact strongly have a strong charge called 
color. 

The rest of the particles in the Standard Model which are not force carriers, 
the fermions, are subdivided into different groups depending on their properties 
(See table IT"^ , reminiscent of the way chemists organized the elements into 



^ Bosons are defined as particles with integer spin in quantum mechanics. 
^ Fermions are defined as particles with half odd fractional spin. 
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the periodic table during the second half of the nineteenth century. Fermions are 
divided into quarks (generically represented by the symbol q) and leptons. The 
main difference between these two groups is that quarks interact through the strong 
force while leptons do not. There are six types of leptons and six types of quarks 
(also called the six quark flavors) which are grouped into three generations. Each 
generation consist of two quarks and two leptons. All three generations replicate 
the same set of force charges, the main difference between generations is the mass 
of the particles (for example, the ratio of the masses of e : : r is 1 : 200 : 3500). 
For each fermion there is an anti-fermion with equal mass and spin and opposite 
charge. 

Table 1-2: Fermion symbols classified into quarks, leptons, and the three genera- 
tions along with the generational common charges. 





First 


Second 


Third 


Electric 


Weak 


Has 




Generation 


Generation 


Generation 


Charge 


Charge 


Color 


Quarks 


u 

d 


c 

,s 


t 

b 


+2/3 

-1/3 


+1/2 

-1/2 


Yes 

Yes 


Leptons 


e 


1^ 


T 



-1 


+ 1/2 
-1/2 


No 
No 



1.1.2 Quantum Chromodynamics 

Quantum chromodynamics (QCD) is the part of the Standard Model that 
describes the strong force. QCD is based on local gauge invariance and color 
symmetry. There are three possible color charges for quarks called r (red), b (blue), 
and g (green). Anti-quarks have opposite colors called f, b, g. The strong force 
between quarks only depends on their colors and is independent of their flavor. 

A very important characteristic of the strong force is that gluons themselves 
carry a color charge and an anti-color charge so they can interact with other quarks 
through the strong force and change their color. Since there are three colors and 
three anti-colors one might think that there are nine different gluons. However, 
there is one linear combination of color anti-color states that has no net color and 



leaves a quark unchanged. There are therefore 9 — 1 = 8 gluons. The 8 individual 
gluon color states can be written as follows, 

1 



1 > 



4 > 



7> 



> 



x/2 



[rb + hr) 



3 > = —1= (rf — bb) 
v2 



V2 

—i 



{rg + gr) 



5 > = -jJ^rg-gr) 
6> = ^{bg + gb) 



{b-g - gb) 



—I 

-j={rf + gg - 2bb) 



[1.1) 



and represent the eight different gluons that exist in nature. The single and unique 
color state left out is called a color singlet, 



|9 >= —={rr + bb + gg) 
v3 



[1.2) 



which is invariant under a redefinition of the color (a rotation in color space). In 
group theory this decomposition of the color states into an octet and a singlet is 
denoted by 3 (8) 3 = 8 © 1. It is worth noting here that a colorless sate, such as |3 > 
or 1 8 >, is not necessarily a color singlet state. 

This situation is analogous to the perhaps more familiar example of two spin 
1/2 particles. Each particle can have their spin up {]) or down (|) along the z axis 
corresponding to four possible combinations reprsented by each giving a total spin 
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= or 1 represented hy \S >■ The S = 1 states form a triplet, 



1 +1 > 



TT> 



1 



(I Ti>+I iT>) 



(1.3) 



1 > 



V2 



1 



1 > 



and there is a singlet state with S 



1 



(in>-iiT>)- 



(1.4) 



> 



V2 



Since gluons themselves carry color they can interact with each other through 
the strong force. This interaction among the force carriers is unique to the strong 
force and, when included in perturbative QCD calculations, leads to two important 
properties observed in nature called "asymptotic freedom" and "color confinement." 
Asymptotic freedom means that the interaction gets weaker at short distances. 
Color confinement is the requirement that observed states have neutral color, or in 
other words, they must be in a color singlet state. 

Confinement explains why free quarks or free gluons, which have a net color 
charge, have never been observed. It also explains why no fractional charged 
particles made from a qq bound state have never been observed since it is not 
possible to construct a color singlet for such a state (in group theory terms 
3 ® 3 = 6 © 3 where we have a sextet and a triplet, but no singlet). On the other 
hand, color singlets can be constructed for a qq or qqq system. The color singlet for 
qq is simply the state shown on the right side of Equation 11.21 The color singlet for 
qqq can be obtained from the decomposition 3(8)3®3 = 10©8©8©1 and is. 



Particles that are bound states of qq are abundant in nature and are called 
mesons, those that are bound states of qqq are called baryons and also abound in 



qqq > color singiet= —r={rgb - grb + brg - bgr + gbr - rbg). 
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nature. Both groups are collectively called hadrons. There have been hundreds of 
different hadrons observed, confirming the validity of QCD and the quark flavors. 
Nevertheless, QCD leaves room for more possibilities. A bound state of two gluons 
gg (sometimes called a glueball) can be in a color singlet and in principle could be 
observed. 

Although glueballs are allowed by QCD there is no convincing experimental 
observation of one. Another possibility are bound states that are a mixture of the 
previous states, such as aqq + 13 gg with arbitrary a, /?. These states are called 
"hybrid mesons" or simply "hybrids". It is believed jl], [21, [3] that hybrids are 
necessary to explain the spectrum of light mesons between 1-2 GeV. 

As with any quantum mechanical theory, the spectrum of bound states is 
a fundamental test. ^ Glueballs are allowed by QCD, yet there is no conclu- 
sive experimental evidence of their observation, despite intense experimental 
searches jH], [Hj complemented by lattice QCD calculations 0, jH! and other 
theoretical contributions like bag models [9 , fiux-tube models ^U], QCD sum 
rules weakly bound bound-state models [T2j, and QCD factorization formalism 
models [13]. Physicists cannot be sure they understand QCD until such states are 
observed, or until they can explain why we cannot observe them. 
1.1.3 Introduction to the Radiative Decays of the T(1S) 

The T(1S) is a meson composed of hh quarks. Mesons of this kind, composed 
of a quark and an anti-quark of the same flavor, are in general called quarkonia. 

There is a convention behind the name of the T(1S) . The T symbol is 
reserved for particles composed of hb where the combined spin of the quark and 
anti-quark is 1. The "IS" symbol is borrowed from atomic spectroscopy with the 



^ A good example of this is the successful description of the hydrogen atom's 
spectrum by quantum mechanics. 
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"1" meaning that the bb pair are in the lowest-energy bound-state, and the "S" 
meaning that the bb have a relative angular momentum L = 0. A description of 
particle naming conventions can be found in |14j . 

The T(1S) is unstable, existing for only about 10~^^ seconds after which it 
decays into daughter particles (which in turn decay themselves if they are not 
stable). The term "radiative decay" is reserved to any T(1S) decay where one of 
the stable daughters is a photon. 

The different ways the T(1S) can decay must obey the symmetries in nature. 
Examples of such symmetries are "parity" and "charge conjugation", both of which 
will be described soon. Symmetries are very important in the standard model. The 
usefulness of symmetries can be seen in Noethers' theorem, which states that for 
every symmetry there is an associated conserved quantity. Examples of Noethers' 
theorem are the conservation of momentum connected to translational invariance 
and the conservation of angular momentum associated with rotational invariance. 

In the next sections we will use the parity and charge conjugation symmetries 
to find T(1S) decays which conserve the associated symmetry constants and are 
allowed by nature. In particular, we will show that the radiative decay of the 
T(1S) through a photon and two gluons (see Figure IT"^ is allowed. The key 
observation is that the two gluons must be in a color singlet since both the T(1S) 
and the radiated photon have no color and color must be conserved. This means 
that the two gluons satisfy color confinement and could form a glueball, although 
more conventional meson states, or hadrons in no bound state at all, are also 
possible outcomes. Regardless of what the gluons do, their energy will eventually 
manifests itself as hadrons. Sometimes, two charged hadrons of opposite charge will 
emerge. This work examines those two hadrons from a radiative T(1S) decay to 
experimentally probe the gg spectrum. 
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Pcirity 

The parity operator, P, reverses the sign of an object's spatial coordinates. 
Consider a particle \a > with a wave function ^a(a;, t). By the definition of the 
parity operator, 

P^{x,t)=Pa^a{-X,t) (1.6) 

where Pa is a constant phase factor. If we consider an eigenfunction of momentum 

^p(x,t) ^e'^-^-^'^ (1.7) 

then 

P^p-(f , t) = Pa'^pi-X, t) = Pa^-pix, t), (1.8) 

SO that any particle at rest, with p — 0, remains unchanged up to a multiplicative 
number. Pa, under the parity operator. States with this property are called 
eigenstates with eigenvalue Pa- Pa is also called the intrinsic parity of particle a, 
or more usually just the parity of particle a, with the words at rest left implicit. 
Since two successive parity transformations leave the system unchanged, = 1, 
implying that the possible values for the parity eigenvalue are Pa — ±1. 

In addition to a particle at rest, a particle with definite orbital angular 
momentum is also en eigenstate of parity. The wave function for such a particle in 
spherical coordinates is, 

^r.lm{S,t)^Rr,l{r)Yr{e,<P), (1.9) 

where (r, 9, (j)) are spherical polar coordinates, Rni{r) is a function of the radial 
variable r only, and the Yi^{9, (f)) is a spherical harmonic. 

The spherical harmonics are well known functions which have the following 
property, 

vne, 0) = i-i^n-K -e,7r + <!>). (1.10) 
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Hence 

P^nlm{x,t) = Pa^nlm{-X,t) = P„(-l)'^,;^(f, (1.11) 

proving that a particle with a definite orbital angular momentum I is indeed an 
eigenstate of the parity operator with eigenvalue Pa(— 1)'- 

The parities of the fundamental fermions cannot be measured or derived. All 
that nature requires is that the parity of a fermion be opposite to that of an anti- 
fermion. As a matter of convention fermions are assigned P = +1 and anti-fermions 
are assigned P — —1. In contrast, the parities of the photon and gluon can be 
derived by applying P to the field equations resulting in = — 1 and Pg — —1. 

The T(1S) has P = PbPii-l)^ = -1 since L = 0. 

Parity is a good quantum number because it is a symmetry of the strong and 
electromagnetic force. This means, that in any reaction involving these forces, 
parity must be conserved. 
Chcirge Conjugation 

Charge conjugation is simply the operation which replaces all particles by 
their anti-particles. In quantum mechanics the charge conjugation operator is 
represented by C. For any particle |a > we can write, 

C\a>^Ca\a> (1.12) 

where is a phase factor. If we let the C operator act twice to recover the original 
state |a >, 

\a >— C'^\a >— C{ca\a >) = CaC\a >— CaCa\a > (1-13) 

which shows that CaCa — 1. If (and only if) a is its own anti-particle, it is an 
eigenstate of C. The possible eigenvalues are limited to C = Ca = Cg, = ±1. 

All systems composed of a the same fermion and an anti-fermion pair are 
eigenstates of C with eigenvalue C = (— 1)^^+'^). This factor can be understood 
because of the need to exchange both particles' position and spin to recover the 
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original state after the charge conjugation operator is apphed. Exchanging the 
particles' position gives a factor of (—1)^ as was shown in the previous section, 
exchanging the particles spin gives a factor of (— l)*^"^^ as can be verified by 
inspecting Equations 11.31 and II .41 and a factor of (-1) which arises in quantum field 
theory whenever fermions and anti-fermions are interchanged. With this result we 
can calculate the charge conjugation eigenvalue for the T(1S) and obtain C = —1 
since L + S = 1. 

The photon is an eigenstate of C since it is its own anti-particle. The C 
eigenvalue for the photon can be derived by inserting C into the field equations and 
is Cj = —1. 

Finally, we consider a system composed of two gluons that are in a color 

singlet. The two gluons are bosons and they must have a symmetric wave function, 

under a gi ^ g2 exchange. Under this exchange, the orbital angular momentum 

part of the wave function contributes with a factor (—1)'^, the spin part of the 

wave function for two spin 1 particles contributes with a factor of (—1)"^, and the 

color singlet part of the wave function contributes with a factor of +1 since it is 

symmetric. To ensure that is symmetric we need L + S* to be even. This implies 

that C = (—1)'^+'^ = +1. The L = and L = 1 possible gg bound states are shown 

in Table The J^'" = 1^+ is peculiar because it is impossible for a qq system to 

have these quantum numbers. If this state is ever observed, it must be a glueball. 

Experimental searches for such a state have been done, for example, in fSj- 

Table 1-3: Possible gg bound states with L = or L = 1. The possible quantum 
numbers are limited by the condition that J = L + S he even, which is needed to 
ensure a symmetric wave function for the two gauge bosons. 



L 


s 










0++ 





2 


2++ 


1 


1 


0-+,l-+,2-+ 
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Charge conjugation is a symmetry of the strong and electromagnetic force. For 
those particles that are eigenstates of (7, C is a good quantum number because in 
any reaction involving these forces C must be conserved. 

The P and C values of various particles used in this analysis are shown in 
Table II 41 along with their quark composition, orbital angular momentum, and 
internal spin. 

Table 1-4: Symbol, name, quark composition, angular momentum (L), internal spin 
(S), parity (P), and charge conjugation eigenvalues (C) for a few of the particles 
used in this analysis. 



Symbol 


Name 


Quark Composition 


L 


S 


P 


c 


T(1S) 


Upsilon(lS) 


bb 





1 


-1 


-1 


7r+ 


Pion 


ud 








-1 


X 




Pion 


du 








-1 


X 


K+ 


Kaon 


us 








-1 


X 


K- 


Kaon 


su 








-1 


X 


P 


Proton 


uud 





1/2 


+1 


X 


P 


Anti-proton 


uud 





1/2 


-1 


X 


7 


Photon 


X 


X 


1 


-1 


-1 



Possible Decays of the T(1S) 

At this point we can understand the different possible ways the T(1S) can 
decay. The Possible decays are limited because the strong and electromagnetic 
force must conserve color, C and P. 

The simplest possibility is for the bb pair to interact electro-magnetically and 
annihilate into one virtual photon ^ . This is allowed by parity, charge conjugation 
symmetry and color conservation. The decay of the T(1S) to one gluon is not 
allowed by color conservation and is therefore forbidden. T(1S) decays to two 
photons are forbidden by charge conjugation. T(1S) decays to two gluons are 



^ Such a photon is called virtual because it cannot conserve the 4-momentum of 
hb and is unstable, only existing for a brief period of time, as allowed by the uncer- 
tainty principle, after which it decays. 
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also forbidden by charge conjugation. T(1S) decays to 3 gluons are allowed (3 
gluons can form a color singlet). T(1S) decays to 3 photons are also allowed, but 
are largely suppressed by decays to one photon since 3 successive electromagnetic 
interactions are much less likely to occur than a single one. Finally, T(1S) decays 
two one photon and two gluons are also allowed under the condition that the two 
gluons be in a color singlet state. 

The three different possible T(1S) decays with least amount of interactions 
(also called lowest order decays) are shown in Figure IT^ 




Figure 1-1: Lowest order decays of the T(1S) allowed by color conservation, charge 
conjugation symmetry, and parity, (a) Shows the decay into three gluons, (b) shows 
a radiative decay, and (c) shows the electromagnetic decay through a virtual pho- 
ton that in turn decays electromagnetically into a pair of charged fundamental 
particles, such as quarks or charged leptons (the charged leptons are represented by 
the symbol /). 

Observable Resonances 

In this work we search for a resonance X produced in a radiative T(1S) decay. 



T(l^) ^ 7X 



(1.14) 
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Using charge conjugation symmetry on both sides, 

-1 = -ICx. (1.15) 

Therefore, Cx = +1- 

In order for us to observe X it must decay into two charged hadrons, 

X^h+h-, (1.16) 

where h = 7T,K,p are the hadrons whose momentum we are going to measure. 
Applying charge conjugation to this last decay, 

+1 = (-1)^+^ (1.17) 

where L and S are respectively the angular momentum and spin of the h'^h~ 
system. This last equation implies that L + S must be even. This has consequences 
for the possible X parities we can observe. By parity conservation in ll.l(i[ 

P^ = (-l)^(-l)(-l) = (-l)^ (1.18) 

For h = IT, K, S = and L = J must be even, which implies that Px = +1. For 
h = p, S can be or 1 and particles with both positive and negative parities can be 
detected. 

Table 1-5: Possible S, L, J, P and C values for X from the radiative decay 
T(IS') 7X reconstructed in different decay modes. 



Decay Mode 


S 


L 


J 


P 


C 


X TT + TT- 





even 


even 


+1 


+1 


X K+K- 





even 


even 


+1 


+1 


X —>■ pp 





even 


even 


+ 1 


+1 


X ^ pp 


1 


odd 


even and odd 


-1 


+1 
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1.2 Radiative Decays of Quarkonia Overview 

Theoretical models exist for glueball production in quarkonia decay jT^l 
and for the glueball spectrum. For example, a quenched lattice calculation [7] 
predicts a J^^ = 2++ glueball ^ in the 2.2 GeV/c^ mass region (see Figure IT"^ . 
According to Table II 51 in the charged pion and kaon modes, we are limited to 
detect glueballs in the leftmost column where P = C = +1 of Figure H"^ while in 
the proton mode we are restricted to the two left most columns. 

12 
10 
8 

E 6 

o 

4 
2 








++ 



-+ +- 

PC 



Figure 1-2: Quenched lattice calculation result for the glueball spectrum for differ- 
ent P and C values. The mass scale is shown in terms of a scaling parameter from 
the QCD lattice calculation named vq on the left and in absolute terms on the right 
by taking r^^ = 410 MeV. 



^ Here J stands for the internal angular momentum (spin) of the glueball, 
J = L + S. 
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For some time bound states not considered to be pure glueballs, such as 
/2(1270), /4(2050), 7], and i]' have been observed in J/ip radiative decays at 
the 10"'^ production level jTHj ^ . In 1996, the BES collaboration claimed the 
observation of a resonance they called the /j(2220) particle in the radiative J/ip 
system at the 10~^ level ^7] (see Figure IT"^ . A lot of excitement was generated at 
the time because it is possible to interpret the /j(2220) as a glueball. However, this 
result has not been confirmed. 




2.2 2.4 
Mass (GeV) 



Figure 1-3: The mass spectrum obtained by the BES collaboration in radiative 
J/ip decays into different hadronic modes. 



Key to identifying a particle as a glueball are (a) suppressed production in 
two-photon collisions (unlike quarks, gluons don't carry electric charge and do not 



^ The J/ip particle is a cc bound state. Since the strong force is flavor blind the 
situation in J/ip decays is in principle similar to that in T(1S) decays. 
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couple to photons), and (b) flavor symmetric decays, since a pure glueball has no 
valence quarks. 

CLEO has already done several studies of radiative decays of the T(IS'). 
Naively, one expects these types of decays to be suppressed by a factor 

[(gb/gc)K/m;,)]2 ^ 0.025 (1.19) 

with respect to J/ ip radiative decays. This comes from noticing that the quark- 
photon coupling is proportional to the electric charge and the quark propagator 
is approximately 1/m for low momentum quarks. In 1999, CLEO made the first 
observation of a two-body T(IS') radiative decay JH]- The spin of the observed 
resonance could not be measured, but its mass and width where consistent with the 
/2(1270) particle. Under this assumption, comparing the measured branching ratio 
of T(1S) 7/2(2220) to the measured branching fraction of the J /ip — > 7/2(1270), 
a suppression factor of 0.06 ± 0.03 was obtained. After the BES result for the 
/j(2220) in radiative J/ip decays, a corresponding search was performed by CLEO 
in the radiative T(IS') system ^Hl- This analysis put limits on the /j(2220) 
production in radiative T(1S) decays. 

In this work we are privileged to have available the largest collection of 
radiative T(1S) decays in the world. With it, we can study the structure of color 
singlet gg hadronization, and shed more light on the /j(2220) result from BES. 



CHAPTER 2 
EXPERIMENTAL APPARATUS 

To carry out our study of the di-hadron spectrum we need to first produce 
the T(1S) resonance and secondly observe its daughter particles fiying away at 
relativistic speeds. These two tasks are respectively accomplished by the Cornell 
Electron Storage Ring and the CLEO III detector. 

2.1 The Cornell Electron Storage Ring 

The Cornell Electron Storage Ring (CESR), located at Cornell University, is a 
circular particle accelerator that produces e~^e~ collisions. 

In order to produce such collisions electrons and positrons need to be created, 
accelerated and stored. CESR's different components, shown in Figure have 
been carrying out this task since 1979. 

A typical CESR run begins at the linear accelerator (LINAC) where electrons 
and positrons are produced. To create positrons, electrons are evaporated off a 
filament and linearly accelerated by electromagnetic fields towards a tungsten 
target. The collision creates a spray of electrons, positrons and photons. The 
electrons are cleared away with magnetic fields and the positrons are introduced 
into the synchrotron. The filling procedure is identical, except that the tungsten 
target is removed. 

Once the electron and positron beams are introduced into the synchrotron, 
they are accelerated to the operating energy. In hour case they are accelerated to 
the point where their combined energy is the T(1S) mass, 9.46 GeV. 

Once the beams are at the desired energy they are transfered to the storage 
ring, where they will remain for about an our. At one point in the storage ring the 
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beams are forced to cross paths. This is the point where e+e coUisions occur 
and where the center of the CLEO III detector is located. If the accelerator is 

N 



Figure 2-1: The Wilson Laboratory accelerator located about 40 feet beneath 
Cornell University's Alumni Fields. 

performing well a high collision rate results. A high collision rate is crucial for the 
success of an accelerator and the experiments it serves. The important figure is the 
number of possible collisions per second per unit area; this is called the luminosity. 
In order to maximize the luminosity, the beams are focused as small as possible at 



^ This point is not fixed in space, but varies from event to event inside a small 
volume of space called the interaction region (IR). 
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the IR. During the CLEO III installation several magnetic quadrupoles were added 
to CESR to improve the beam focus. CESR has consistently outdone itself while 
collecting luminosity over the years (see Figure . 




Figure 2-2: CESR yearly luminosity. The gaps in 1995 and 1999 correspond to 
down times when the CLEO II. V and CLEO III detectors where being installed. 

The CLEO III detector measures the time integrated luminosity over a period 
of time by counting how many times a benchmark process occurs during a certain 
time interval at the IR. For redundancy, there are two benchmark processes that 
are used, one where the e^e~ particles interact to produce a new e^e— pair, 
and the other one where the e^e~ annihilate and produce two photons. Using 
the known cross-section for each process, the number of events is converted to a 
luminosity called the Bhabha integrated luminosity for the first process, and the 77 



20 

integrated luminosity for the second one. The term "integrated" is sometimes left 
out and the total luminosity is referred to as simply the Bhabha or 77 luminosity, 
with the time integration left implicit. 

2.2 The CLEO III Detector 

When the e^e~ collision occurs, the two particles are annihilated we enter 
the world of particle physics. Nature decides what to do with the energy from 
the annihilation. We have no chance of directly observing what is happening at 
the annihilation point, but eventually long lived semi-stable particles are created 
that fly off at relativistic speeds. These particles carry information about what 
happened after the e+e~ collision, and can tell us what nature did. The CLEO III 
detector has the important mission of detecting and measuring such particles. 

As one can see in Figure l2-3[ the CLEO III detector is a composite of many 
detector elements. These sub-detectors are typically arranged as concentric cylin- 
ders. The entire detector is approximately cube shaped, with one side measuring 
about 6 meters, and weighs over 500 thousand kilograms. 




Magnet Barrel Muon 

Ifon Chambers 



Figure 2-3: The CLEO III detector. 
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As implied by the name, CLEO III is not the only CLEO detector. CLEO 
III was preceded by CLEO II.V, CLEO II, CLEO I.V and the original CLEO 
detectors. The CLEO III detector was a major upgrade compared with the 
previous version of CLEO [201 , and has an improved particle identification 
system together with a new drift chamber and a new silicon vertex detector. 

2.2.1 Superconducting Coil 

All the detector subsystems except for the muon chambers are located inside 
a superconduction coil. The coil remains unchanged since CLEO II. It is kept in 
a superconducting state by liquid helium. The purpose of the coil is to provide a 
1.5 Tesla magnetic field, which is uniform to 0.2%, to bend the paths of charged 
particles in the detector. By measuring how much much a charged particle bends, 
experimenters can measure the momentum of the particle. 

The coil inner radius is 1.45m and its outer radius is 1.55 m, with a radial 
thickness of 0.10 m. The total length of the coil in z is 3.50 m. It is wound from 
a 5mm x 16 mm superconducting cable (Al surrounding Cu-NbTi strands). It 
is wound in 2 layers, with 650 turns per layer, on an aluminum shell. When in 
operation a current of 3300 amps fiows through the coil. 

2.2.2 Tracking System 

After particles from the interaction point pass through the beam pipe, they 
begin to encounter the active detector elements of the tracking system. There are 
two sub-detectors responsible for tracking the curving path of charged particles. 
The first one encountered by particles is the silicon vertex detector, and the second 
one is the central drift chamber. The CLEO III tracking system is responsible 
for tracking a charged particle's path and measuring its momentum. Typical 
momentum resolution is 0.3% (1%) for 1 GeV (5 GeV) tracks. The tracking system 
also measures ionization energy losses with an accuracy of about 6%. 
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Silicon Vertex Detector 

The silicon vertex detector in CLEO III [221; also called SVD III, is a four- 
layer barrel-only structure with no endcaps that surrounds the beam pipe. This 
detector (see Figure 1^"^ provides four cf) and four z measurements covering over 
over 93% of the solid angle. The average radius of inner surface of the four layers is 
25 mm, 37.5 mm, 72 mm, 102 mm. The detector is constructed from 447 identical 
double-sided silicon wafers, each 27.0 mm in 0, 52.6 mm in z and 0.3 mm thick. 
The wafers are instrumented and read out on both sides. The instrumentation on 
each side consists of an array of aluminum strips on the wafer surface. These strips 
are connected to preamplifiers at the end of the detector. The inner side has 512 
strips in the z direction and the outer side has 512 in the direction. Therefore 
each wafer contains 512+512 sensors. The 447 wafers are arranged in the 4 layers, 
as follows: 7 sections in 0, each with 3 wafers in z, total = 21 wafers in the first 
layer; 10 sections in 0, each with 4 wafers in total = 40 wafers in the second 
layer; 18 sections in 0, each with 7 wafers in z, total =126 wafers in the third 
layer; 26 sections in phi, each with 10 in Z, total = 260 wafers in the fourth layer. 

Charged particles traversing the wafer lose energy and create electron hole 
pairs. Approximately 3.6 eV is required to create a single electron-hole pair. The 
electrons and holes then travel in opposite directions in the electric field applied 
to the surfaces of the wafers until they end up on the aluminum strips, and the 
detector registers a "hit" . When combined together, the hit on the inner side of a 
wafer and the hit on the outer side give a measurement of the (z, 0). The wafer 
position itself determines r. 

The Central Drift Chamber 

The CLEO 111 central drift chamber (DR III) is full of a gas mixture with 60% 
Helium and 40% propane held at about 270 K and at a pressure slightly above one 
atmosphere. The drift chamber is strung with array of anode wires of gold-plated 
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Figure 2-4: View of the SVD III along the beampipe. 

tungsten of 20 fim in diameter and cathode wires of gold-plated aluminum tubes 
of 130 /im in diameter. The anodes are kept at a positive voltage (about 2000 V), 
and the cathodes are kept grounded, which provides an electric field between the 
anode and the cathode wires. Anode and cathode wires are often called "sense" 
and "field" wires respectively. 

As a charged particle passes through the DR III, it interacts electromagneti- 
cally with the gas molecules giving energy to the outer electrons which become free 
in a process called ionization. The free electrons from the ionized gas molecules 
drift in the electric field toward the nearest anode wire. As the electrons get close 
to the anode, the electric field becomes very strong which causes an avalanche 
as further ionization is induced. The result of the avalanche is a large number of 
electrons collapsing upon the sense wire in a very short amount of time (less than 
one nanosecond) . When this happens to a sense wire, we say that there is a "hit" . 
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The current on the anode wire from the avalanche is amphfied and collected at 
the end of the anode wire. Both the amount of charge and the time it takes it to 
move to the end of the detector are measured. A calibration of the drift chamber 
is used to convert the amount of charge to a specific ionization measurement of the 
incident particle. A calibrated drift chamber can also convert the time to roughly 
measure the position along the sense wire where the charge was deposited. 

The wires are strung along the z direction. About 2/3 of the outer part of 
the drift chamber (the farthest part from the interaction point) is strung in with a 
slight angle (about 25 miliradians) with respect to the z direction to help with the 
z measurement. Wires strung in the z direction are called "axial" wires, while those 
that are strung at a slight angle are called "stereo" wires. 

The DR III consists of an inner stepped section with 16 axial layers, and an 
outer part with conical endplates and 31 small angle stereo layers. There are 3 
field wires per sense wire and they approximately form a 1.4 cm side square. The 
drift resolution is around 150 yum in r — and about 6 mm in z. All wires are held 
at sufficient tension to have only a 50 /im gravitational sag at the center (2; = 0). 
There are 1696 axial sense wires and 8100 stereo sense wires, a 9796 total. 
2.2.3 Ring Imaging Cherenkov Detector 

The Ring Imaging Cherenkov (RICH) detector [23] is a new detector subsys- 
tem for CLEO III. It replaces the CLEO II. V time of flight system designed to 
measure particles' velocities. 

Cherenkov radiation occurs when a particle travels faster than the speed of 
light in a certain medium, 

V > c/n. (2.1) 

Where v is the velocity of the particle, c is the speed of light in vacuum, and n 
is the index of refraction of the medium the particle is traveling in. The charged 
particle polarizes the molecules of the medium, which then turn back rapidly 
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to their ground state, emitting radiation. The emitted hght forms a coherent 
wavefront if f > c/n and Cherenkov hght is emitted under a constant Cherenkov 
angle, 5, with the particle trajectory forming a cone of light. The cone half-angle is 
given by the Cherenkov angle which is, 

cos (5= — = -^. (2.2) 

vn pn 

If the radiation angle, S, is measured, the speed of the incident particle is known. 
This measurement, combined with the momentum measurement from the tracking 
system, gives a measurement of the particles mass, and can be used in particle 
identification. 

The threshold velocity at which Cherenkov radiation is emitted is Vmm = 
When a particle traveling at the threshold velocity transverses the medium a very 
small cone with 5 ~ is produced. The maximum emission angle occurs when 
Vmax = c and is given by 

cos6^ax = -- (2.3) 
n 

The RICH (see Figure consists of 30 modules in phi, 0.192 m wide and 
2.5 m long. The detector starts at a radius of 0.80 m and extends to 0.90 m. Each 
module has 14 panes of solid crystal LiF radiator at approximately 0.82 m radius, 
0.192 m wide, 0.17 m long, 1 cm thick. Inner separation between radiators is 
typically 50/im. The LiF index of refraction is n = 1.5. The radiators closest to z 
= in each module have a 45 degree sawtooth outer face, to reduce total internal 
refiection of the Cherenkov light for normal incident particles (see Figure 1^"^ . The 
radiators are followed by a 15.7 cm (radial) drift space filled with pure N2. The 
drift space is followed by the photodetector, a thin-gap multiwire photosensitive 
proportional chamber. 

With this index, particles in the LiF radiator with beta = 1 produce 
Cherenkov cones of half-angle cos~^(l/n) = 0.84 radians. With a 16 cm drift 
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Figure 2-5: The RICH detector subsystem. 

space, this produces a circle of radius 13 cm. The RICH is capable of measuring 
the Cherenkov angle with a resolution of a few miliradians (see Figure ITTj) . This 
great resolution allows for good separation between pions and kaons up to about 3 
GeV as Figure 12-81 shows. 
2.2.4 Crystal Calorimeter 

The CLEO Crystal Calorimeter (CC) is composed of 7784 thallium-doped 
Csl crystals. Each crystal is 30 cm long (16.2 radiation lengths) with 5cm x 5cm 
square front face. The crystals absorb any incoming electron or photon which 
cascades into a series of electromagnetic showers. The electronic system composed 
of 4 photo-diodes present at the back of each crystal are calibrated to measure the 
energy deposited by the incoming particle. Other incoming particles other than 
photons and electrons are partially, and sometimes fully, absorbed by the crystal 
giving an energy reading. 

The CC is arranged into a barrel section and two endcaps, together covering 
95% of the solid angle. The CC barrel section is unchanged since CLEO II; the 
endcaps have been rebuilt for CLEO HI to accommodate the new CESR interaction 
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Figure 2-6: The two kinds of RICH LiF radiators. For normal incidence particles 
(2; ~ 0) a sawtooth radiator is necessary to avoid internal reflection. 

region quadrupoles. The barrel detector consists of an array of 6144 crystals, 128 
in (f) and 48 in z, arranged in an almost-projective barrel geometry. That is, the 
crystals are tilted in z to point to a few cm away from the interaction point, and 
there is also a small tilt in 0. The CC barrel inner radius is 1.02 m, outer radius 
is 1.32 m, and the length in z at the inner radius is 3.26 m. It covers the polar 
angle range from 32 to 148 degrees. The barrel crystals are tapered towards the 
front face (there are 24 slightly different tapered shapes), the endcap crystals are 
rectangular, but shaved near the outer radius to fit in the container. The CC 
endcaps consist of two identical end plugs, each containing 820 crystals of square 
cross-section, aligned parallel to the beam line (not projective). There are 60 
crystals in the "fixed" portion of the "keystone" piece of the endcap. 760 in the 
part that slides. The keystone is made up of two parts, one on top that has 12 
crystals that for mechanical removal reasons is separate from a container holding 48 
crystals. The endcap extends from 0.434 m to 0.958 m in r. The front faces are z = 
±1.308 m from the interaction point. It covers the polar angle region from 18 to 34 
degrees in +z, and 146 to 162 in -z. 
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Figure 2-7: Cherenkov angle resolutions per track as a function of radiator row for 
Bhabha events. Row 1 corresponds to the two rows closest to 2; = 0, etc. 

The photon energy resolution in the barrel (endcap) is 1.5% (2.5%) for 5 
GeVphotons, and goes down to 3.8% (5.0%) for 0.1 GeVphotons. 
2.2.5 Muon Detectors 

The muon detectors (MU) are the most external subsystem of the CLEO III 
detector. They remain unchanged from CLEO II, and are composed of plastic 
proportional tubes embedded in the magnet iron return yoke. They cover 85% 
of the An solid angle (roughly 30-150 degrees in polar angle). If a series of hits is 
detected in the muon chamber layers they most likely correspond to muons because 
other particles are blocked by the iron. Besides detecting muons, the heavy iron 
of the return yoke protects the inner subsystems of the CLEO III detector from 
cosmic ray background (except for cosmic ray muons of course). 
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Figure 2-8: Pion fake rate as a function of particle momentum for kaon efficiency of 
80% (circles), 85% (squares) and 90% (triangles). 

There are three planes of chambers in the barrel section, arranged in 8 octants 
in (j). The plastic barrel planes lie at depths of 36, 72, and 108 cm of iron (at 
normal incidence), corresponding to roughly 3, 5, and 7 hadronic interaction 
lengths (16.8 cm in iron) referred to as DPTHMU. There is one plane of chambers 
in each of the two endcap regions, arranged in 4 rough quadrants in cj). They lie 
at z =± 2.7 m, roughly covering the region 0.80 < | cos(6')| < 0.85. The planar 
tracking chambers use plastic proportional counters at about 2500 V with drift gas 
of 60% He, 40% propane, identical to (and supplied by the same system as) the 
drift chamber gas. Individual counters are 5 m long and 8.3 m wide, with a space 
resolution (along the wire, using charge division) of 2.4 cm. The tracking chambers 
are made of extruded plastic, 8cm wide by 1 cm thick by 5 m long, containing eight 
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tubes, coated on 3 sides with graphite to form a cathode, with 50 /im silver-plated 
Cu-Be anode wires held at 2500 V. The orthogonal coordinate is provided by 
8cm copper strips running perpendicular to the tubes on the side not covered by 
graphite. 



CHAPTER 3 
ANALYSIS OF THE DATA 

This work builds on the techniques developed by previous CLEO radiative 
T(1S) analyses [HI, [El, modified to be used with CLEO HI data. A new 
technique based on kinematic fitting is developed to, together with the new 
RICH detector, improve efficiency and particle identification. 

We search for radiative T(1S) decays in the modes T(1S) '^tt^tt^ , 
and 'ypp. The e+e^ collision data has both resonant events, where the e~^e~ 
annihilate to give a T(1S) , and continuum events, where the e^e~ collision does 
not give a T(1S) . To be sure we are observing T(1S) and not a continuum process, 
the continuum must be subtracted by using a pure source of correctly scaled 
continuum events. Pure continuum data can be obtained by operating CESR at an 
energy different form the T(1S) mass. After subtracting the underlying continuum, 
we examine the di-hadron invariant mass spectrum in search of resonances. We 
determine the spin and production helicity (the projection of the spin on the 
momentum vector at production time) of any found resonances by examining the 
photon and hadron angular distributions. 

3.1 Data Sample 

The analysis presented here is based on CLEO HI data. Throughout this 
document, and unless otherwise stated we use 77 luminosity (the luminosity types 
used by CLEO are defined at the end of Section 2.1). We prefer to use the 77 over 
the Bhabha luminosity because the resonant process T(1S) —>■ e^e~ artificially 
increases the reported Bhabha luminosity by about 3% in T(1S) data. This extra 
contribution would need to be accounted for when doing a continuum subtraction. 
Choosing the 77 luminosity avoids this complication. 

31 
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The CLEO III data is divided into numbered sets. Sets 18 and 19 have a 
luminosity of 1.13 ± 0.02/6^^ in the beam energy range 4.727-4.734 GeV. This 
data, which we call the T(1S) data (or simply the IS data), has both resonance 
e+e^ T(1S) and continuum events. We take the number of resonant events 
from |24j, Nr{is) = (2.1 ± 0.1) x 10''. This number contrasts with the previous 
generation measurement of CLEO II, where A^t(is) ~ 0.15 x 10'' were available. 

In datasets and 18 and 19 there are also 0.192 ± 0.004/6"^ taken below the 
T(1S) beam energy (4.714-4.724 GeV). This data, which we call the T(lS)-off data 
(or simply the IS-off data), has relatively low statistics and corresponds to purely 
continuum events. 

To improve our continuum statistics we use 3.49 ± 0.07/6^^ from datasets 
9, 10, 12, 13, and 14 of data taken near the center-of-mass energy of the T(4S), 
which for our purposes is defined as data with beam energy in the 5.270-5.300 GeV 
range. This set of data, which we call the T(4S) data (or simply the 4S data), is a 
source of pure continuum because no T(4S) —>■ BB resonant event can survive our 
"cuts" ^ (the cuts are presented in Section 3.3). 
3.1.1 Continuum Subtracted Distributions 

We use the continuum data taken at the T(4S) energy to subtract the under- 
lying continuum present in the T(1S) data. This is important because continuum 
background processes like e^e^ 7p with p —>■ tt+tt^, e+e^ —>■ 70 with 
(j) K~^K~ , and direct e+e^ '-/h^h^ (we will use the convention h = tt, K, p 
from now on), look like the signal events we are searching for. To first order, the 
cross section of these continuum process scales like 1/s, where s is the square of 



^ Cuts are simply conditions that an event must satisfy to be considered in the 
analysis. Cuts are necessary to eliminate background that would otherwise make a 
measurement difficult or even impossible. 
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the center of mass energy of the e~^e~ system. Taking the luminosity-weighted 
average beam-energies of each interval, and the 77 luminosities (see Table I3-1|) . we 
calculate that the T(4S) data scaled down by a factor of 0.404 ^ represents the un- 
derlying continuum in the T(1S) data. This is true up to differences in momentum 
distributions and phase space. The error in the continuum scale factor is unknown 
because the luminosity ratio is expected to have a small but undetermined sys- 
tematic error. We make the somewhat arbitrary decision to retain three significant 
digits in the continuum scale factor because it is sufficient for our purposes and 
there are 0.5% effects from second order terms in the cross section formulae. 

To eliminate the contribution of continuum events from a T(1S) data variable 
distribution (e.g., the invariant mass of two tracks, the photon angular distribution) 
we proceed as follows, 

• Obtain the T(4S) data distribution for the same variable. 

• Efficiency correct both the T(1S) and T(4S) distributions, using a 
GEANT [25J based Monte Carlo (MC) simulation of the detector. Examples of 
MC efficiency distributions are shown in Figures f3-3l and 17-31 



^ If we are to be mathematically strict we should calculate the scale factor as 

-^IS run I Sy^^ l~-AS run 

ie run ' .c. "^45 run 

lb runs 4A runs 

with obvious notation. This is equivalent to redefining the average energy as 

1 1 ^ ^ ^run 

runs ' 

However, our energy intervals are sufficiently narrow and this calculation does not 
change the last significant digit of the original scale factor. Similarly, taking into 
account second order terms in the energy dependence of the cross section, like the 
one that appears as m^/ s in the explicit formula for the e^e^ 7p cross section, 
also has an insignificant effect on the scale factor. 
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• Subtract the T(4S) distribution from the T(1S) distribution using a T(4S) 
scale factor of 0.404. 

We call this set of steps "continuum subtraction" by definition. For com- 
pactness, we call such a distribution the "continuum subtracted <variable> 
distribution/plot" or "<variable> continuum subtracted distribution/plot". Except 
for statistical fluctuations and phase space effects, the resulting IS distribution 
should not have any contribution from continuum processes that scale as 1/s. 

It is important to notice that any continuum subtracted distribution is 
efficiency corrected. This means that a fit to a continuum subtracted distribution 
(for example, the continuum subtracted invariant mass distribution) gives the 
efficiency corrected number of events directly. Strictly speaking, this number of 
events is only correct if all the other variables we cut on have the same initial 
distribution in data and MC, or if the efficiency does not depend on them. In this 
note we use "fiat MC" , defined as MC that is generated with a fiat distribution in 
the mass and the helicity angles 6^, 9,, (these angles are defined in the appendix). 

For example, the number of events obtained from the continuum subtracted 
invariant mass fit needs to be corrected to account for the fact that the helicity 
angle distributions are not flat in data (see for example Figure l4-9t and the 
efficiency is highly dependent on these variables (see Figure 17"^ . This correction is 
done in Section 7.2. 

A summary of the results from this section is shown in Table 1311 

3.2 Event Selection 

Event selection for T(IS') jh'^h~ is straightforward and can be thought of 
in terms of three major stages. 

First we skim the data, keeping only those events that contain exactly one 
high-energy photon and two tracks. Next, we require that the total 4-momentum 
of these three elements be consistent with the colliding e'^e~ 4-momentum. Finally, 
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Table 3-1: Summary of the data used in this analysis. The continuum scale factor 
is obtained using 77 luminosities because the Bhabha luminosity is artificially high 
during T(1S) running due to the process T(1S) — > e'^e~ . 



T(1S) T(4S) T(lS)-off 

Dataset 18, 19 9, 10, 12, 13, 14 18, 19 

Average Eb^am{GcV) 4.730 5.286 4.717 

Range of E6eam(GeV) 4.727-4.734 5.270 - 5.300 4.714-4.724 

£(e+e-)(/6-i) 1.20 ±0.02 3.56 ±0.07 0.201 ± 0.004 

>C(77)(/&"^) 1.13 ±0.02 3.49 ±0.07 0.192 ± 0.004 

T(1S) continuum scale factor 1 0.404 5.84 



in the third stage, we project the surviving data onto the three different hadronic 
modes via hadron separation and QED suppression cuts. 

However, checking the 4-momentum involves using the tracks masses. This 
means that the information from stage 2 should somehow be useful in stage 3. This 
is indeed the case, and the details of how we do it are revealed in this section. 
3.2.1 Skim Cuts 

We skim over the data in the "hardGam" subcoUection. The hardGam 
subcoUection was developed with this type of analysis in mind. For an event that 
passes the triggers ^ to be classified as hardGam it must pass the following cuts, 

• eGaml > 0.5 

• eSh2 < 0.7 

• eOverPl < 0.85 

• eVis > 0.4 

• aCosTh < 0.95 where "eGaml" is the highest isolated shower energy 
relative to the beam energy, "eSh2" is the energy of the second highest shower 
relative to the beam energy, "eOverPl" is the matched shower energy relative to 



^ Triggers are basic criteria that an event must satisfy to recorded during the 
data collection processes. Triggers are designed to get rid of trash and noise and 
reduce the size of the data sample while keeping all of the important information. 
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the momentum of the track with highest momentum, "eVis" is the total energy 
detected (charged tracks are assumed to be pions) relative to the center of mass 
energy, and "aCosTh" is the z component of the unitary total momentum vector. 

Monte-Carlo predicts that about 75% of the generated T(15') •yh'^h^ signal 
passes the hardware and software triggers and gets classified as hardGam. 

As mentioned above, we use this data to make our skim. To write an event 
from the hardGam subcollection into our skim we require the following topological 
cuts, 

• There are exactly two "good tracks" ; there can be any number of tracks that 
are not "good tracks" but these are not used in the analysis. We define a "good 
track" as a track that satisfies the following cuts; drift chamber track ionization 
energy loss (dE/dX) information is available, the ratio of number of wire hits to 
those expected is between 0.5 and 1.5, the pion fit has x^/d-o.f. < 20 (here d.o.f. 
stands for degrees of freedom), and the distance of closest approach to the beam 
spot in the x-y plane (called DBCD) is less than 5 — 3.8P (mm) if P < 1 GeV/c 
(where P is the tracks momentum in GeV/c) and less than 1.2 mm for tracks with 
P > 1 GeV/c. This DBCD cut is common in the more sophisticated CLEO II/II.5 
analyses. It performs better than a simple DBCD < 5 mm cut, because it takes 
into account the fact that traks with higher momentum have a better measurement 
of DBCD since they scatter less. 

• There is exactly one "good shower", there can be any number of showers 
that are not "good showers" but these are not used in the analysis. We define a 
"good shower" as an unmatched shower with energy > 4 GeV. 

These topological cuts are about 85% efficient for generated signal events that 
have passed the triggers and have been classified as hardGam. 

The overall skim efficiency is between 60-65%, depending on the mode (see 
Table l3~4|l . 
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3.2.2 Analysis Cuts 

After our skim we call any cuts we make "analysis cuts" . These cuts are done 
at analysis time and are mode dependent. As a convention, and unless otherwise 
stated, efficiencies for individual analysis cuts are reported relative to the events in 
the skim (not relative to the events generated). 

4-momentum Cut 

All fully reconstructed events should have the 4-momentum of the e^e~ 
system. This constraint is usually implemented with a simple two-dimensional 
AE-p box cut, where AE is the difference between the reconstructed energy for 
the event and the colliding e^e~ energy (Ecm), and p is the magnitude of the 
reconstructed total momentum for the event. Typical values for these cuts are 
-0.03 < AE/EcM < 0.02 and p < 150 MeV/c (taken from fT^). 

The traditional AE-p box cut is somewhat useful. However, it does not take 
into account the correlation between the measured energy and momentum. Indeed, 
the signal lies in diagonal bands in the AE-p plane, making a box-shaped cut not 
optimal (see Figures 13-lb and 13-lb ). 

We use an alternative approach to the 4-momentum cut. After a simple 
substitution, = = Ecu — E^+h- (where E^ is the photon's energy, p^ is 
the magnitude of the photons momentum, Eqm is the energy of the e^e~ system, 
and -E/i+zi- is the energy of the hadron pair), we can write the E-p conservation 
equations as: 

Ph+h- + {EcM - Eh+h-)p^ = PcM (3.1) 

where, Ph+h- is the di-hadron momentum, pj is the photon's momentum unit 
vector, and pcM is the momentum of the e"*"e~ system (which is a few MeV because 
of the crossing angle). Equation 13. II is a 3-constraint subset of the 4-momentum 
constraint and has the convenient property of avoiding the use of the measured 
photon's energy, which has non-Gaussian asymmetric errors. It is important to 
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notice that Equation 13.11 contains the di-hadron energy, therefore, it can help 
discriminate between the various particle hypotheses. 

We proceed as follows. After vertexing the hadron pair using the beam spot 
with its error matrix, we calculate the photon's direction from the hadron pair's 
vertex and the shower position. We then fit the event to the 3 constraints expressed 
in Equation (1) using the techniques outlined in and cut on the of the 
3-constraint fit, Xe-pW < ^OO- '^^ complete the 4-momentum requirement, we 
calculate AE{h) = E^+h- + — Ecm, where Eh+h- is the updated di-hadron 
energy after the constraint, and E^ is the measured photon's energy, and require 
—0.050 * EcA4 < AE{h) < 0.025 * Eqm- Furthermore, we now have available 
Xe-pW differences between different particle hypotheses, which help in particle 
identification (ID). This is discussed in more detail in the next section. Figure 
compares the performance of the old and new approaches to the E-p cut. 

At this point it is a good idea to check that the 4-momentum cut rejects 
background events that make it through our skim cuts (see Section 3.1). These 
events typically have one high energy shower, two tracks, and (an) additional 
element(s). One such background is T(IS') — > ■yn^n^n^. Out of 25000 T(IS') 777, 
with 1] tt+tt^tt'^ MC events ^ only 4 survive our 4-momentum cut. We conclude 
that our 4-momentum cut is good at rejecting background events that pass the 
skim cuts but have additional elements such as an extra photon, 7r°, pair of tracks, 
etc. 



^ Thanks to Vijay Potlia for generating these events 
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Figure 3-1: Distributions for different 4-momentum cuts. Signal MC events 
(T(IS') 77r+7r^) are represented by the black dots, and "background" MC events 
(e^e^ — > 70, — K^K^) appear as (red) triangles. Plot a) has no cuts. Plot b) 
has the old 4-momentum box-cut. Plot c) has the new 4-momentum cuts. Plot d) 
has the new 4-momentum cuts, and also a cut defined as xl;-p(^) ~ XE-pi^) < 0- 
The particle ID potential of the newly available Xe-pW evident. 

Hadron Separation 

We define Ax|^(/ii — for the particle hypotheses hi and h2 of our charged 
track pair (e.g. Axj^(7r — K)) as follows, 

-2\og{CRicH{ht)) + 2\og{CRicH{ht)) (3-2) 
-2\og{CRicH{K)) + 2\og{C Ricuih^)). 

where the idea is to combine the dE/dX and RICH information into one number. 
Pairs of tracks with Axj^jihi — /i2) < are more likely to be of type hi than of type 
h,. 
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In practice, we only add a tracks RICH information if its momentum is above 
the Cherenkov radiation threshold for both mass hypotheses and there are at least 
3 photons within 3a of the Cherenkov angle for at least one of the mass hypothesis. 
We also require that both hypotheses were actually analyzed by RICH during 
pass2. 

In addition to RICH and dE/dX, and as hinted in the previous section, the 
difference in Xe-pW from the constraint expressed in Equation 13. II can help the 
particle ID (see Figures 13 Ik and 13-ll i). We define, 

Axl.,{h - h,) = xl-p{hi) - Xl-p{h2). (3.3) 

Events with Ax^_p(/ii — /i2) < are more likely to be of type hi than of type /i2. 

In this analysis, to select hi and reject /i2 the default cut is simply Ax|^(/ii — 
/i2) < 0. This simple cut is highly efficient, has low fake rates, and is sometimes 
sufficient. However, out of the six possible cases when one hadron fakes another, 
there are three important cases where it pays off to also use Ax^^pihi — h2) 
together with Axjoi^i ~ ^2) in an optimal way; 

1. vr background to K. This background comes from the continuum process 
e^e^ 7P, p TT+TT^. 

2. TT background to p. Again, this background comes from the continuum 
process e+e~ 7p, p tt+tt^. 

3. K background to p. This background comes from the continuum process 
e^e~ — s> 70, — * K^K^ . 

In other words, the important cases occur when the lighter mass hypothesis 
fakes the heavier mass hypothesis. 

Mathematically, one would expect that simply adding both Axjo and Axe-p 
together (like we just did when combining RICH and dE/dX), and cutting on the 
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grand Ax^ is the way to go. Unfortunately, this simple approach fails because of 
large non-mathematical tails in the individual Xe-p distributions. 

Instead, for each of these three cases we define the best cut values (ci, C2) in 
AxjD < Cl and Axe-p < C2 as those that maximize 



F{Ci,C2) 



(3.4) 



^Re + W*Rf 

where Re (-R/) is the efficiency (fake rate) of the particle ID cuts and W is the 
rough ratio of the background to signal in the data sample for each case. Each W 
can't be known a priori, but a rough idea of its value can be obtained by doing 
a first iteration of the analysis with, for example, Ci = C2 = 0. We use W = 20, 
W = 60, and W = 30 for cases 1-3 respectively. Figure I2E21 shows F(ci,C2) for each 
case. Table 1^21 shows the optimized cut values and their effect on particle ID. 




Figure 3-2: Contour plots of F(ci,C2) as defined in Equation 13.41 for different sig- 
nals and backgrounds from flat MC. Top has a K signal and a vr background. Cen- 
ter has a p signal and a vr background. Bottom has a p signal and a K background. 



Table 3-2: Cut values, efficiencies, and fake rates for Axjd < Ci and Ax|;_p < C2 in flat MC. The different Ci, C2 values are 
chosen so F from Equation 13.41 is maximun. Efficiency and fake rates of each cut are reported relative to events in the skim. 
Statistical errors in the efficiencies are 0.1% or less. Errors in the fake rates are statisticall and are shown for completeness 
only. 







Cut Value 


Signal efficiency (%) 




Fake rate (%) 








Cl 




Axjo AxLp Both 






Both 


TT faking 


K 


3 


-2 


94.2 95.2 89.9 


1.73 ±0.09 


14.5 ±0.3 


0.31 ±0.04 


TT faking 


P 


-1 


-2 


99.0 99.8 98.9 


0.89 ±0.06 


1.03 ±0.07 


0.03 ±0.01 


K faking 


IP 


-1 


3 


98.8 99.3 98.1 


1.89 ±0.09 


1.09 ±0.07 


0.10 ±0.02 
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Quantum Electrodynamic Suppression 

Quantum Electrodynamic (QED) background in this analysis comes from the 
abundant processes e+e^ — > 7e+e^ and e+e~ — ^ 7yU+yU^. 

To reject e+e^ 7e+e^ we require both tracks to have a matched shower with 
energy E such that \E/p{'n-) — 0.95| > 0.1, and to have Axj^ih — e) < 0. 

To reject e+e^ — * 'yfi^fi^ we simply require DPTHMU < 5 [DPTHMU 
was defined in Section 2.2.5) for both tracks in the K and p modes because particle 
ID cuts make the pion (and therefore the muon) fake rate small (see Tables 
and l3-fi|l . For the vr mode we cannot use particle ID in a practical way because 
muons and pions have similar masses. Instead, to separate pions and muons we 
use a much stronger cut requiring that both tracks be within the barrel part of 
the muon chamber (| cos(0)| < 0.7), both have P > 1 GeV/c and both have 
DPTHMU < 5^ . 

To improve the overall muon suppression cut efficiency with virtually no 
increase in muon fakes, we fiag an event as "not muonic" if any of the tracks 
deposit more than 600 MeV in the CC. This increases the cut efficiency by about 
90% in the vr mode and makes the detector more hermetic. 
3.2.3 Cut Summary, Efficiencies, and Fake Rates 

Table 13-31 summarizes the cuts used in this analysis. Figure 13 31 shows the 
overall Monte Carlo efficiency after all cuts. Figure l3^ shows the fake rates 
according to the MC and the data for different particle ID cuts. The data fake 
rates and limit fake rates, which are measured for pions and kaons faking other 



5 Other analyses (for example [IH]) typically use DPTHMU < 3, our CLEO III 
MC has too many pion tracks with 3 < DPTHMU < 5 and two few tracks with 
DPTHMU = 0. 
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hadrons, are calculated from the p and peaks in the continuum (see Figures 14-21 
and I5-1|) . Tables 13 41 through 13-61 summarize the results of Figures 13 31 and 13 41 



Table 3-3: Cuts used in this analysis. 



Motivation 


Definition for T(15') —>■ '-/h^h (where h = tt,K, or p) 


Skim cuts 


Data acquisition 


Event must pass hardware (Levell) and software triggers 
(pass2) and be of type hardGam 


Topological 


There are only two good tracks and only one isolated shower 
with E>A GeV 


Analysis cuts 


Reconstructed event must 
have 4-momentum of the 
center-of-mass system 


Xl-p{h) < 100.0 and -0.050 < AE{h)/EcM < 0.025 


Hadron separation 


Default is, Axj^ < 0. The three cases with a large fake pop- 
ulation because of continuum processes use a simultaneously 
optimized cut on ^Xid ^^"^ ^Xe-p^ ^"^^ summarized in 
Table lO 


QED background 


Both tracks have a matched shower energy E that satisfies 
\E/p{tt) - 0.95| > 0.1, and Ax?^(/i - e) < 


QED background 
e+e^ 7/1"*" /i~ 


For /i = K, p both tracks have L'PTiJMf/ < 5. For 
h = 7T (At least one track has a matched shower energy > 
600 MeV) or (( Both tracks have cos(e) < 0.7 and P > 
1 GeV/c) and (both tracks have DPTHMU < 5)) 
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Table 3-4: Efficiencies in % for cuts (as outlined in Table for flat signal MC. 
Efficiencies in the second group are reported relative to the number of candidates 
that make it to the skim. The third part of the table shows the overall reconstruc- 
tion efficiency. Statisticall errors are 0.1% or less. 



Skim cuts 


TT 


K 


P 


Hardware Trigger 


89.7 


89.0 


90.9 


Software Trigger 


96.5 


96.2 


96.6 


hardGam 


78.4 


79.9 


76.5 


Topological {ritracks = 2 Sz > AGeV) 


73.3 


67.6 


72.3 


Overall Skim efficiency 


64.3 


60.3 


62.4 


Analysis cuts 


TT 


K 


P 


4-momentum 


98.6 


98.5 


99.0 


QED e+e^ 7e+e^ suppresion 


93.9 


87.4 


93.1 


QED 7yU+/i^ suppresion 


74.7 


93.0 


98.3 


Hadron separation 


97.1 


89.0 


97.8 


Overall analysis efficiency 


66.9 


79.1 


89.1 


Overall reconstruction efficiency 


43.0 


47.6 


55.6 



Table 3-5: Efficiencies in % for cuts (as outlined in Table for flat 4S MC. Effi- 
ciencies in the second group are reported relative to the number of candidates that 
make it to the skim. The third part of the table shows the overall reconstruction 
efficiency. Statisticall errors are 0.1% or less. 



Skim cuts 


TT 


K 


P 


Hardware Trigger 


88.4 


88.0 


89.0 


Software Trigger 


94.6 


94.6 


94.8 


hardGam 


79.2 


80.5 


77.9 


Topological {utracks = 2 Sz > AGeV) 


75.2 


69.5 


75.8 


Overall Skim efficiency 


66.9 


62.7 


66.6 


Analysis cuts 


TT 


K 


P 


4-momentum 


97.1 


97.0 


97.4 


QED e+e^ ■ye'^e' suppresion 


93.4 


89.1 


93.2 


QED — > 7yU+/i~ suppresion 


75.8 


92.5 


98.1 


Hadron separation 


98.0 


91.4 


96.6 


Overall analysis efficiency 


67.4 


80.2 


87.3 


Overall reconstruction efficiency 


45.0 


50.3 


58.2 
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Table 3-6: Final efficiencies and fake rates after all cuts in %. Statistical errors 
for efficiencies are 0.1% or less. Statistical errors for fake MC rates are shown for 
completeness only. The IS DATA corresponds to the IS off resonance data sample. 
The MC was generated flat. 





TT cuts 


K cuts 


p cuts 


IS TT MC 


43.0 


0.14 ±0.01 


< 0.007 


4S TT MC 


45.0 


0.29 ±0.02 


< 0.01 


IS p DATA 




< 1.2 


< 0.2 


4S p DATA 




0.20 ±0.06 


< 0.06 


IS K MC 


1.27 ±0.03 


47.6 


< 0.02 


4S K MC 


0.92 ±0.04 


50.3 


< 0.06 


IS (j) DATA 


< 3.8 




< 2.0 


4S (/) DATA 


4.14 ±0.69 




< 0.4 


IS p MC 


0.08 ±0.02 


0.34 ± 0.04 


55.6 


4S p MC 


0.05 ±0.01 


0.37 ±0.01 


58.2 
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Figure 3-3: Final efficiencies for each mode as a function of invariant mass for the 
IS (solid) and 4S (dashed) Monte Carlo data. The MC was generated with a flat 
angular distribution. 
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Figure 3-4: Fake rates for IS (hollow circles) and 4S (hollow squares) according 
flat MC. The (red) downward pointing triangle is obtained using the p peak in 
data. 



CHAPTER 4 

EXCLUSIVE RADIATIVE DECAY T{1S) 771+71" 
In Figure the 7r"''7r~ invariant mass plot is shown for both IS and 4S data. 
Figure shows the continuum subtracted 7r+7r^ invariant mass distribution (as 
defined in Section 3.2.1) with the most likely statistical fit overlayed (which is 
described in the next section). The number of events within IF of the p region 
(0.62 - 0.92 GeV/c^) left after the continuum subtraction is 200 ± 300, and 50 of 
these belong to the /o(980) low- mass tail. 

4.1 Robustness of The Mass Distribution 
In j2Z] the decay T(1S) — > 77r°7r° is analyzed, and it is shown how the analysis 
stream warps the shape of the reconstructed resonance. This effect, which arises 
because of the particular 7r° behavior, raises problems when fitting the invariant 
mass distribution. 

In Section 3.2.1 we claimed that if the data and MC had the same 6.y and 
distributions, the fit to the continuum subtracted invariant mass distribution 
automatically gives the correct efficiency corrected number of events. 

Here we test this claim. To this end, we generate 10000 T(1S) 7/2(1270), 
with /2(1270) 7r+7r^ with fiat 6.y and distributions. We treat this MC 
as data and carry out the first two steps of the continuum subtraction process. 
The resulting peak has a mass of 1.278 ± 0.002 GeV/c^ and a width of 0.193 ± 
0.006 GeV/c^, consistent with the generated mass and width of 1.275 GeV/c^ and 
0.185 GeV/c^. More importantly, the number of reconstructed events from the fit is 
10040 ± 180, which is consistent with the number of generated events. 

Figure 14-11 shows the reconstructed and efficiency corrected events, a fit to 
them, and the generated events. 
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We conclude that there is no warping of the mass distribution, and that the 
analysis stream behaves like we expect when obtaining the efficiency corrected 
number of events from data. 



■ Reconstructed ond Eff. Corrected 
I 1 MC Generated 

g^o"! 




1.00 1.10 1.20 1.30 1.40 1.50 

7t*7t" Invariant Mass 

Figure 4-1: Reconstructed events and efficiency corrected events, a ffi to them, and 
the generated events for T(1S) 7/2(1270), with /2(1270) tt+tt" with flat 9^ 
and distributions. 

4.2 Statistical Fit of the Invariant Mass Distribution 

The results of this section are summarized in Table FTl 
Figure shows possible signals for for T(IS') 7/o(980), T(IS') — >• 
7/2(1270) and T(IS') 7/4(2050). To determine the number of events in each 
signal we ffi the invariant tt+tt^ mass continuum subtracted distribution with three 
spin-dependent, relativistic Breit-Wigner line shapes. The masses and widths are 
allowed to float, except for the width of the /4(2050) which has a very large error 
if allowed to float and is set to its PDG value. The PDG values ^H] for the mass 
and width of the three resonances are mj(j(98o) = 980 ± 10 MeV/c^, rjQ(98o) = 
70 ± 15 MeV/c2, mf^^i270) = 1275 A ± 1.2 MeV/c^, Tf^^i27o) = 185. 1^^^ MeV/c^, 
^/4(2050) = 20 34 ± 11 MeV/c2, and r/,(2050) = 222 ± 19 MeV/c^. 

To measure the statistical significance of each signal we do multiple fits 
fixing the signal area to different values while letting the mass and width of 



51 

the signal whose significance is being measured to float within 2a of the PDG 
values. At the same time the rest of the flt parameters are flxed to the values that 
originally minimized the x^. We assign each of these multiple flts a probability 
proportional to e~-^^/^ and then normalize. We calculate the chance of the signal 
being due to a random fluctuation by adding the normalized probabilities for the 
flts with a negative or signal. This method fails for the highly signiflcant /2(1270) 
signal because the e~-^^/^ value of flts with negative or /2(1270) signal is of the 
order of 10~^^ and our computing software can only handle numbers as small as 
10~^^. For completeness we state that the signiflcance of this signal is < 10~^^. 

To measure the upper limit for the /j(2220) we also do multiple fits 
for different flxed signal values, while keeping its mass and width constant at 
^/j(2220) = 2.234 GeV/c2 and r/j(222o) = 17 MeV/c^ as in [HI . The resulting 
probability plot is shown in Figure FHl 

We flnd clear evidence for the /2(1270), evidence for the /o(980) and weak 
evidence for the /4(2050). We also put a 90% confldence level upper limit on 
/j(2220) production. Fit results are shown in Table FTl 

Table 4-1: Results for T(IS') 77r^7r^. The branching fractions of /2(1270) and 
/4(2050) are taken from the PDG JB]- Errors shown are satistical only. 



Mode 


Area 


B.F. (10-^) 


Signiflcance 


7/o(980),/o(980) -^TT+TT- 




-■-•0-0.6 


8.3 X 10-6 (4.3a) 


7/2(1270) 


1230 ± 100 


10.4 ±0.8 


< 10-"^^ (> 14cr) 


7/4(2050) 


85 ±30 


3.6 ± 1.3 


5.2 X 10-3 (2.6(t) 


7/j(2220),/j(2220)^7r+7r- 


< 13 


< 6.2 X 10-2 





4.3 Angular Distribution of The Signal 

The helicity angle distributions of 9^^ and ^7 are deflned and described in the 
appendix. In this section we flrst obtain the helicity angular distributions of the 
/o(980), /2(980), and /4(2050) and then flt them to the predictions of the helicity 
formalism. 
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In practical terms, obtaining the helicity angle distribution of a particular 
resonance from data consists of two steps. First, we choose an invariant mass 
interval around the resonance peak to select events from the resonance and obtain 
a helicity angular distribution which has both signal and background events. 
Second, we subtract the contribution to the helicity angular distribution of the 
background events in the chosen mass interval to obtain what we want; the helicity 
angle distribution of the resonance. 

Choosing the mass interval is not a trivial thing. If its too wide there will 
be too much background, and if its to narrow there will not be enough signal. 
To choose the optimum mass interval we need to know how much signal and 
background we are selecting. Therefore, the two steps described in the previous 
paragraph are related. How we deal with this is revealed in the next two sections. 
4.3.1 Optimum Mass Interval 

As described above, the first step in obtaining the helicity angle distribution 
for a resonance is choosing an invariant mass interval to select events from such a 
resonance. A standard IF (which corresponds to 1.6a for a spin Breit-Wigner) 
cut around the mean mass of the resonance can be chosen as a "standard" interval. 
We could proceed this way, but in our case because of the large subtractions 
involved when obtaining the angular distribution, a considerable increase statistical 
significance of each bin in the helicity angular distribution can be achieved by 
choosing the mass interval carefully (see the last column of Table . 

Let's consider for example the /o(980) angular distribution. We begin with 
the T(IS') and the T{4S) efficiency corrected distributions. Before the continuum 
subtraction each bin in the /o(980) angular distribution has contributions from 
the high end mass tail of the p and the low mass tail of the /2(1270). In order to 
get the final angular distribution, both of these contributions are taken away by 
first doing a continuum subtraction using the scale factor a = 0.404, and then 
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by subtracting the /2(1270) distribution outside the mass interval scaled by an 
appropriate scale factor f3 equal to the ratio of /2(1270) inside the mass interval 
and outside of it ^ . After all this, the contributions to the bins in the final angular 
distribution are, 

A^/o(980) = nfo(gsO) + ^p(T{15)) + %2(1270) - a^p(T(4S)) " /5^/2(1270) ~ '^^Mm (4-1) 

where ^/o(980)(^/2(i27o)) is the number of /o(980)(/2(1270)) obtained by integrating 
the fitted spin-dependent, relativistic Breit-Wigner function inside the mass 
interval, np(x{l5))(^^p(T(4S'))) is the number of p's inside the mass interval from 
continuum events at the T(1S')(T(4S')) energy, and ^^2(1270) is the number of 
/2(1270) outside the mass interval being used to subtract the contribution of 
/2(1270) to the /o(980) angular distribution. Following this last definition, P = 

'^/a (1270) /'^/a (1270)- 

Each of these terms has an associated error. Assuming that An = ^/n, that 
the efficiency correction has infinite statistics, and ignoring the errors on the 
continuum scale factor a, the overall error is. 



AA^/o(980) = \Jnfo(980) + (a + a^)^p(T{45)) + (2 + + /92)n/2(i27o) 

We arrive to the conclusion that the mass interval (m — Am, rh + Am) which 
produces the helicity angle distribution with smallest relative bin-errors is the one 



^ Strictly speaking this is correct only to first order. The /2(1270) distribution 
itself has a small contamination from /o(980) and /4(2050). This effect, which we 
call cross-contamination, is ignored in Equation 14.11 Later, in section 4.2.2, we will 
show how to eliminate this small cross-contamination using all the resonances. The 
difference in how we actually get the helicity distributions and Equation 14. II has an 
insignificant effect when calculating the optimal mass interval 
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that maximizes, 

V^/o(980) + (a + a^)^p(T(45)) + (2 + P + P^)n/2(i270) 

The plot of FyQ(98o)(?Ti, Am) is shown in figure FUl 

This same technique can be apphed to the /2(1270) and /4(2050) resonances. 
Results are shown in Table 14-21 



Table 4-2: Mean masses, widths in GeV/c^ and inverse of the average bin relative error (F) from background subtractions 
for the angular distribution of different resonances. Standard mean masses and widths, corresponding to IF, are taken from 
the fit in Figure and are labeled with the subscript "s" , while those that maximize F are labeled with the subscript "m" . 
The last column shows the factor by which the effective statistics increase. 



Resonance 


{ms,Ams) 


F{ms,Ams) 


{rhm, Anim) 


F{mm,Amm) 


F'''{mm,A.'mm) 
F^{m,s, Arris) 


/o(980) 


(0.970,0.070) 


5.7 


(0.985,0.060) 


6.0 


1.11 


/2(1270) 


(1.270,0.120) 


18.3 


(1.590,0.420) 


20.6 


1.27 


/4(2050) 


(2.120,0.220) 


3.2 


(2.240,0.250) 


3.4 


1.13 
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4.3.2 Background Subtraction 

The continuum subtraction of the helicity angle distribution is defined in 
Section 2.1. The subtraction of the tails from other resonances requires a closer 
look. Let us call the continuum subtracted helicity angle distribution of the events 
in the /o(980), /2(1270), and /4(2050) mass intervals c/o(98o), c/2(i27o), and 0/4(2050) 
respectively. Let us call the helicity angle distribution of the events that come 
exclusively from the resonance we are trying to select in the same mass intervals 
(that is, the true helicity angle distribution of the resonance) t/(,(98o), t/2(i27o); and 
^/4(2050) respectively. The continuum subtracted helicity angle distribution of a 
resonance is being contaminated by the tails of the other resonances. Keeping this 
in mind we write ^ 

C/o(980) = i'^/oCgSO) + /5^/2{1270) 
^ C/2(1270) = ^/o(1270) + 7^/o(980) + '^'^/4(2050) (^'"^^ 
C/4(2050) = ^/o{2050) + ei/2(1270) 

V 

Where the small numbers /?, 7, 5, and e are the ratios of the number of events 
from a resonance in the mass interval where the contamination is taking place to 
the number of events from the same resonance in the mass interval used to select 
it. Using the (mm, Amm) values in Table and the fit in Figure FHl we obtain 
/3 = 8.5 X 10-^ 7 = 6.0 X 10-2, 5 = 9.7 X 10"^, and e = 0.13. 

To obtain the background subtracted helicity angular distributions we simply 
invert the system of equations expressed in Equation 14.31 The solution can be 



2 Here the cross contamination between the /o(980) and the /4(2050) is ignored. 
There is no mathematical problem in including these contamination terms, but the 
system of equations would be more complicated than it needs to be since this type 
of cross contamination is negligible. 
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conveniently expressed as, 
( 

^/2(1270) = i_^!^_5,(c/o(1270) - 7C/o(980) - '^C/4(2050)) 
* i^/oCgSO) = C/o(980) - /5^/2{1270) (^-^ 
^/4(2050) = C/o{2050) - e^/2(1270) 

V 

As a check, if we ignore the second order terms the previous solution becomes, 

i/2(1270) = C/o(l270) - 7C/o(980) " ^C/4(2050) 
^ ^/o(980) = C/o(980) - /?C/2(i270) 
1^/4(2050) = C/o(2050) ~ eC/2(i270) 

V 

which is indeed is the solution when cross-contamination is ignored. 
4.3.3 Statistical Fit of the Helicity Angular Distributions 

For each resonance, we fit the data to the simultaneous cos and cos Q-,, 
helicity angle distribution obtained from data using Equation 14.41 to the helicity 
formalism prediction. Equations EBIEZI projected on each angle and folded in 
opposite directions around their symmetry axis in order to show both distributions 
on the same plot. The corresponding J value of the best fit for each resonance, 
shown in Figures 14-81 through I4-1()| is defined as the J assignment (Ja). We obtain 
Ja = 1 for the /o(980) which is inconsistent with the known spin of the /o(980) 
which is J = 0. For the /2(1270) and the /4(2050) we obtain Ja = 2 and Ja = 4 
respectively, which is consistent with their known spins. 

To have an idea of how well the angular distribution determines J among the 
hypotheses J = 0, 1, 2, 3, 4, we do a statistical fit for each hypothesis and assign 
each one a probability proportional to e"*^^^"'"'^'"''^'-*/^ where d.o.f. are the degrees 
of freedom in the fit. The resulting normalized probability distributions give an 
idea of the assigned J significance and are shown in Figure E-l II In particular. 
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this figure shows that the Ja — 1 for the /o(980) inconsistency can be due to a 
statistical fluctuation. 



File: 7bohr/user1/luis/cleo3/ntp/data_50.ntp 

ID IDB Symb Date/Time Area Mean R.M.S. 

1 1 1 040807/0140 9112. 0.8080 0.2895 

1 2 1 040807/0140 2.1195E+04 0.7647 0.2072 

98765 981 -1 040807/0203 569.1 1.227 0.7650 
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Figure 4-2: Invariant mass of tt'^tt" for IS (top) and 4S (bottom) data. For the 4S 
data we show a blow-up of the mass region 1.5 — 3 GeV/c^ where the p* can be 
seen. 
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Figure 4-3: Continuum subtracted (as defined in Section 3.2.1) invariant mass 
7r"'"7r~ from T(15') — > 777"^ 7r~. 
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Figure 4-4: Blow up of the invariant mass of tt+tt in the /4(2050) mass region, 
including the upper limit on the /j(2220) shown as a dashed hue. 
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File: /tmp_mnt/bohr/userl/luis/cleo3/mn_fit/u1sgammahh/temp.tmp 

ID IDB Symb Date/Time Area Mean R.M.S. 

1 29 1 040806/1845 1.000 6.294 5.008 

1 25 1 040806/1845 0.9001 5.132 3.650 
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Figure 4-5: Normalized probability distribution for different /j(2220) 
nal areas. The shaded area spans 90% of the probability. 
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Figure 4-6: Contour plot of the inverse of the average relative bin-error from back- 
ground subtractions in the /o(980) (top), /2(1270) (middle), and /4(2050) (bottom) 
angular distribution. 
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Figure 4-7: Angular distribution for the excess events in the /o(980) mass region. 
The fit corresponds to J = 0. 
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Figure 4-8: Angular distribution for tiie excess events in tiie /o(980) mass region. 
The fit corresponds to J = 1. 
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Figure 4-9: Angular distribution for the excess events in tiie /2(1270) mass region. 
The fit corresponds to J = 2. 
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Figure 4-10: Angular distribution for tlie excess events in tiie /4(2250) mass region. 
The fit corresponds to J = 4. 
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Figure 4-11: probability distribution for resonances in the /o(980) (top), 
/2(1270) (middle), and /4(2050) (bottom) invariant mass region when only the 
hypotheses = 0, 1, 2, 3, 4 are considered. 



CHAPTER 5 

EXCLUSIVE RADIATIVE DECAY T{1S) -fK+R- 
In figure the K^K~ invariant mass plot is shown for both resonance and 
continuum running. This figure also has an inset showing the (JdE/dx{K) for both 
tracks for events in the 1.1 — 3 GeV/c^ region. This inset indicates that most 
of the events in this mass region have indeed two kaons and that the amount of 
p refiection is small. Furthermore, when the vr+vr" invariant mass is plotted for 
these events only 40 out of 700 can be fit under a p peak (this fit is what we use to 
calculate the pion faking kaon rate in Table 13-fjj) . 

Figure 1^21 shows our fit to the K^K~ invariant mass continuum subtracted 
plot as defined is Section 3.2.1 with the most likely statistical fit overlayed (which 
is described in the next section). The number of events near the region (1.01 — 
1.03 GeV/c^) left after the continuum subtraction is 50 ± 70. 

5.1 Statistical Fit of the Invariant Mass Distribution 
The results of this section are summarized in Table 

From the measurement of the previous section we expect a small contribution 
(~ 50 events assuming no interference) from /2(1270) K^K~ . We do find 
some evidence for /2(1270) K^K" events (110 ± 40) in the fit. We also find 
strong evidence for the resonance /2(1525) and weak evidence for the /o(1710) 
resonance. The /2(1270), and /o(1710) are fitted with their widths fixed to their 
PDG values [T^ because they have large errors if allowed to fioat. The rest of 
the resonances parameters are consistent with their PDG values jTH], which are 
m/^(i525) = 1525 ± 5 MeV/c2, 1^.(1525) = 76 ± 10 MeV/c^, m^„(ino) = 1715 ± 
6 MeV/c^, and r/o(i7io) = 125 ± 10 MeV/c^. We also observe an excess of events in 
the 2 — 3 GeV/c^ region which we can't attribute to any known resonances. 
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Significances of tlie signals of the identified resonances in tlie fit are calculated 
as described in the previous section. The T{1S) 7/j(2220), /j(2220) K+R- 
upper limit is also calculated using multiple fits, except that this time the events 
under the /j(2220) are of unknown origin, so we use a first order polynomial 
allowed to float. The significance of the excess of events in the 2 — 3 GeV/c^ 
invariant mass region is calculated assuming a normal distribution; we simply add 
up the number of events in each bin along with its error. Results are shown in 
Table lO 



Table 5-1: Results for T{1S) -fK+K' 



Mode 


Area 


Branching Fraction (10 ^) 


Significance 


7/2(1270) 
7/2(1525) 


110 ±40 


23 ±8 


5.4 X 


10-^(3.3(7) 


360l«° 


Q+0.9 


< 10 


-^^(> 14(7) 


7/o(1710),/o(1710)^ir+ir- 


75 ±30 


0.35 ±0.14 


7.5 X 


10-^(3.2(7) 


'yK+K-{2-3 GcV/c^) 


220 ± 20 


1.03 ±0.12 




8.8(7 


7/j(2220),/^(2220) ^i^+i^- 


< 10 


< 5 X 10-2 
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5.2 Angular Distribution of The Signal 

In this section we adapt the ideas presented in Section 3.4.3 to the T(IS') — > 
jK^K^ situation. 

The derived statistical errors from the signal and background subtractions can 
be used to calculate the mass interval which best represents the helicity angular 
distribution. The inverse of the expected average relative bin error as a function of 
the mass interval is shown in Figure 15 5t and the mass interval that maximizes it 
are tabulated in Table 15 21 



Table 5-2: Mean masses, widths inmass and inverse of the average bin relative error (F) from background subtractions for 
the angular distribution of different resonances. Standard mean masses and widths, corresponding to IF, are taken from the 
fit in Figure 15-21 and are labeled with the subscript "s" , while those that maximize F are labeled with the subscript "m" . 
The last column shows the factor by which the effective statistics increase. 



Resonance 


{rhs, Arris) 


F{ms, Arris) 


{frim, Arrim) 


F{mm,Amm) 


F^(m,„,Amm) 


/2(1270) 


(1.276,0.185) 


3.5 


(1.300,0.100) 


4.1 


1.37 


/^(1525) 


(1.540,0.085) 


9.6 


(1.565,0.100) 


9.8 


1.04 


/o(1710) 


(1.760,0.125) 


3.1 


(1.780,0.095) 


3.3 


1.13 
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The tails from the resonances contribute to the continuum subtracted hehcity 
distributions, 

( 

Cf 2{1270) = tf 2(1270) + Ptf ^(1525) 
^ C/^(1525) = ^f^{1525) + 7^/2(1270) + f^^/oClTlO) ^^'^^ 
C/o(1710) = i^/o(1710) + et/^(1525) 

V 

Where again the small numbers /?, 7, 6, and e are the ratios of the number 
of events from a resonance in the mass interval where the contamination is taking 
place to the number of events from the same resonance in the mass interval used 
to select it. Using the (mm, Amm) values in Table and the fit in Figure we 
obtain f3 = 1.3 x IQ-^, 7 = 0.14, 6 = 0.11, and e = 0.19. 

The background subtracted helicity angular distributions are, 

^/^(1525) = iZ^^(c/^(1525) " 7C/2(1270) - SCf,,(nio)) 
^ ^/2(1270) = C/2(1270) - /5i/^(1525) ^^'^^ 
i/o(1710) = C/o(1710) - e^/^(1525) 

V 

The the best fit for each resonance and the excess of events in the 2 — 3 GeV/c^ 
are shown in Figures f5-6l through 15-91 The best spin assignment for the /2(1270) is 
Ja = 2, for the /2(1525) it is Ja = 2, for the it is /o(1710), and it is = 1 for the 
excess of events in the 2 — 3 GeV/c^ mass region. The = 2 value for the /o(1710) 
is inconsistent with its known spin. Also, examination of the normalized helicity 
amplitudes for the /2(1270) reveals that they are inconsistent with those obtained 
for the /2(1270) in the tt+tt^ mode. 

The assigned J probability distributions are shown in Figure 15-111 They reveal 
that the inconsistencies in the /2(1270) and /o(1710) are not significant and can be 
attributed to the statistical uncertainty. 
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File: */bohr/user1 /Iuis/cleo3/ntp/data_50.ntp 

ID IDB Symb Date/Time Area Mean R.M.S. 

1 1 1 040807/0219 1274. 1.320 0.4283 

1 2 1 040807/0219 2031. 1.156 0.2585 

1 10 1 040807/0235 705.0 0.2246 1.208 

0.1480 1.170 



600 




Mass (GeV/cO 



Figure 5-1: Invariant mass of K^K~ for IS (top) and 4S (bottom) data. For 
the 4S data the inset shows the (JdE/dx{K) for both tracks for events in the 
1.1 — 3 GeV/c^ mass region. This inset indicates that most of the events are con- 
stant with having two kaons. 
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Figure 5-2: Continuum subtracted invariant mass of K^K~ 



76 



File: Generated internally 



ID 


IDB 


Symb 


Date/Time 


9 


3 


-32 


040807/1422 


98765 


981 


-1 


040807/1426 


98765 


982 


-2 


040807/1426 



Area Mean R.M.S. 

794.8 1 .739 0.4861 

308.2 2.086 0.5922 

318.3 2.091 0.5835 



40 



20 



1.80 



2.00 



2.20 

Mass (GeV/c^) 



2.40 



2.60 



Figure 5-3: Blow up of the /j(2220) region, with the 90% CL upper hmit over- 
laid. The mass and width are taken to be m/^(222o) = 2.234 GeV/c^ and 
r/,(2220) = 17 MeV/c2 asin 1^. 
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File: /tmp_mnt/bohr/userl/luis/cleo3/mn_fit/u1sgammahh/temp.tmp 

ID IDB Symb Date/Time Area Mean R.M.S. 

1 29 1 040807/1426 1.000 4.746 3.582 

1 25 1 040807/1426 0.9054 3.971 2.723 




Figure 5-4: Normalized probability distribution for different /j(2220) 
signal areas. The shaded area spans 90% of the probability. 
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Figure 5-5: Contour plot of the inverse of the average relative bin-error from back- 
ground subtractions in the /2(1270) (top), /2(1525) (middle), and /o(1710) (bot- 
tom) angular distributions. 
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MINUITx Fitto Plot 9&82 

UlS^Gamma h h. -ABS(COS(PHOTHETA)) axis 
File: Generated internally 
Plot Area Total/Fit 517.49/517.49 
Func Area Total/Fit 129.01 / 129.01 

■i= 17.8 for 20- Sd.o.f., 



5-AUG-2004 23:16 
Fit Status 3 
E.D.M. 5.475E-05 

C.L.= 40.4% 



Errors 
Function 1 : J2 
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7.95385E-03 
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Figure 5-6: Background subtracted K^K angular distribution in the /2(1270) 
mass region as defined in the text. The fit corresponds to J = 2. 
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MINUITx Fitto Plot 9&72 

UlS^Gamma h h. -ABS(COS(PHOTHETA)) axis 
File: Generated internally 
Plot Area Total/Fit 284.44 / 284.44 

Func Area Total/Fit 367.87/367.87 

■i= 20.9 for 20 - 3 d.o.f., 

Parabolic 



5-AUG-2004 23:16 
Fit Status 3 
E.D.M. 3.771 E-06 

C.L.= 23.3% 
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Figure 5-7: Background subtracted K^K^ angular distribution in the /2(1525) 
mass region as defined in the text. The fit corresponds to J = 2. 



81 



MINUITx Fitto Plot 9&92 

UlS^Gamma h h. -ABS(COS(PHOTHETA)) axis 
File: Generated internally 
Plot Area Total/Fit 129.47/129.47 
Func Area Total/Fit 79.031 / 79.031 

■i= 28.0 for 20 - 3 d.o.f., 



5-AUG-2004 23:16 
Fit Status 3 
E.D.M. 1.721 E-05 

C.L.= 4.5% 



Errors 
Function 1 : J2 
AREA 

* 

50 



39.516 
2.99860E-04 
1.5696 



Parabolic 

± 10.80 
± 0.3847 
± 0.2633 



Minos 
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O.OOOOE+00 + 0.4645 
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Figure 5-8: Background subtracted K^K angular distribution in the /o(1710) 
mass region as defined in the text. The fit corresponds to J = 2. 
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MINUITx'' Fitto Plot 9&102 
U1S-^Gamma h h. -ABS(COS(PHOTHETA)) axis 
File: Generated internally 
Plot Area Total/Fit 



Func Area Total/Fit 

18.9 for 20 - 

Errors 

Function 1: J1 
AREA 185.43 
0.15542 

100 



41 1 .73 / 41 1 .73 
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2 d.o.f., 



5-AUG-2004 23:16 
Fit Status 3 
E.D.M. 4.735E-09 
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Figure 5-9: K^K angular distribution for events within the 2 — 3 GeV/c^ 
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MINUITx Fitto Plot 9&92 

UlS^Gamma h h. -ABS(COS(PHOTHETA)) axis 
File: Generated internally 
Plot Area Total/Fit 129.47/129.47 
Func Area Total/Fit 62.194/62.194 

■i= 33.3 for 20 - 1 d.o.f., 



5-AUG-2004 23:16 
Fit Status 3 
E.D.M. 3.484E-15 

C.L.= 2.2% 



Errors 
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Figure 5-10: Background subtracted K^K^ angular distribution in the /o(1710) 
mass region as defined in the text. The fit corresponds to J = 0. 
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Figure 5-11: Ja probability distribution for resonances in the /2( 1270) (top), 
/2(1525) (middle-top), /o(1710) (middle-bottom) invariant mass region and the 
excess of events in the 2 — 3 GeV/c^ (bottom) region when only the hypotheses Ja 
= 0, 1, 2, 3, 4 are considered. 



CHAPTER 6 
EXCLUSIVE RADIATIVE DECAY T{lS)^-fpp 

In figure IFPTl the pp invariant mass plot is shown for both IS and 4S data, with 
an inset showing that most of the events indeed have a proton and an anti-proton 
and that the p and refiections are small. The enhancement at 3.1 GeV/c^ in 
the 4S plot corresponds to the process e+e~ iJ/i^ with J/ip ^ pp. This 
enhancement is not as pronounced in the IS plot because the IS data only has 
32.4% of the luminosity of the 4S data. The number of events in the J/ip invariant 
mass region after continuum subtraction is 6 ± 3 and is consistent with 0. 

Figure 16-21 shows the continuum subtracted invariant pp mass distribution 
(as defined in Section 3.2.1) with a 90% confidence level upper limit for /j(2220) 
overlaid. A direct fit to the /j(2220) yields 12 ± 5 events. 

There is an excess of events in the continuum subtracted invariant mass plot. 
We measure this excess and the upper limit of e+e^ 7/j(2220), /j(2220) pp 
the same way we measured the excess of events inside 2 GeV/c^ < m{K^ K^) < 
3 GeV/c^ region and the /j(2220) upper limit in Section 5.1. Results are shown in 
Table lO 

The pp angular distribution for the mass range 2 GeV/c^< m{pp) < 3 GeV/c^ 
is shown in Figure IFT^ with the most likely J assignment, Jq = 1, fit overlayed. 
The probability distribution for is shown in Figure 



Table 6-1: Results for T(IS') •ypp. 



Mode 


Area 


B.F. or 90 % U.L. (10-^) 


Significance 


7PP (2-3 GeV/c^) 


85± 18 


0.41 ±0.08 


4.85(T 


7/j(2220),/j(2220)^pp 


< 20 


< 9 X 10-2 
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File: */bohr/user1/luis/cleo3/ntp/data_50.ntp 

ID IDB Symb Date/Time Area Mean R.M.S. 

1 1 1 040807/0303 104.0 2.334 0.3628 

1 2 1 040807/0303 141.0 2.187 0.3419 

1 10 1 040807/0304 139.0 0.2583 1.259 

0.1224 1.119 



1 

" J 


1 ' 

r 

p 


1 

Jl 


' 1 

hn n 




1 

ir 


' 1 ' 1 


1 


— H- 


— 1 — h- 

1 . 


-h 

I 

[L 
1 


^ 6 
3 

i 2 

a -2 

= -6 

- 


' 1 

6 -2 

OdE/c 


1 ' 
1 

2 6 
lx(P) 

n n n r 


Tini ,n 1 





1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 
Mass (GeV/c^) 



Figure 6-1: Invariant mass of pp for IS (top) and 4S (bottom) running. For tiie 4S 
running the inset is consistent with the events having a proton and an anti-proton. 
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File: Generated internally 



ID 


IDB 


Symb 


Date/Time 


9 


3 


-32 


040807/1446 


98765 


981 


2 


040807/1446 


98765 


982 


1 


040807/1446 



Area Mean R.M.S. 

85.43 2.521 0.3093 

72.17 2.685 0.3291 

55.91 2.813 0.2553 




Figure 6-2: Invariant mass of pp. The plot is continuum subtracted and efficiency 
corrected. An overlay with the 90% confidence level upper limit for /j(2220) 
is shown. The mass and width are taken to be m/j(222o) = 2.234 GeV/c^ and 
r/,(2220) = irMeV/c2 as in 117]. 
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Figure 6-3: Normalized probability distribution for different /j(2220) — >■ pp signal 
areas. The shaded area spans 90% of the probabihty. 
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MINUITx'' Fitto Plot 9&22 
UlS^Gamma h h. -ABS(COS(PHOTHETA)) axis 
File: Generated internally 
Plot Area Total/Fit 256.69 / 256.69 

Func Area Total/Fit 1 34.57 / 1 34.57 

19.5 for 20- 2d.o.f., 
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Function 1: J1 
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Figure 6-4: angular distribution for the excess of events in the mass range 
2 GeV/c2< m{vv) < 3 GeV/c^ 
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Figure 6-5: Ja probability distribution for the excess of events in the mass range 
2 GeV/c^< m{pp) < 3 GeV/c^ when only the hypotheses Ja = 0, 1, 2, 3, 4 are 
considered. 



CHAPTER 7 
SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are any sources of experimental uncertainty other 

than the statistical ones. Limits on the accuracy of our detector simulation and any 

physical processes that interferes with the experimental measurement are typical 

examples of systematic uncertainties. These uncertainties need to be identified, 

quantified, and if possible, corrected. 

7.1 Cuts 

From now on we report individual cut efficiencies relative to the events that 
survive all the other cuts. 

According to our MC, most of the skim cuts, except for the hardCam re- 
quirement, are nearly 100% efficient. Therefore, such cuts should not be a source 
of systematic uncertainty. The only skim cut worth taking a closer look at is the 
hardCam cut, which is about 93% efficient in MC. Measuring this efficiency in our 
data is not possible because events not classified as hardGam are not in the data 
to begin with. Closer examination of our MC reveals that that most of the 7% 
inefficiency in hardGam comes from the eOverPl cut (4%), some from the Sh2 cut 
(2%), and the rest (1%) from the other cuts present in hardGam. We can measure 
the efficiency of eOverPl in data by looking at the eOverP2 distribution of data 
tracks from p and decay. We are satisfied by this check on the MC modeling (see 
Tables 17-11 and I7-2|) . and we won't measure how well MC models the rest of the 
hardGam cuts, which are 97% efficient in MC. 
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To quantify the quality in the MC modehng of our analysis cuts in the vr 
and K modes we use the p and signals present in our data^ . In this case the 
process is straightforward. We measure the p and signal signals in data and MC 
over a floating background function with all cuts in place, and with all in place 
cuts except the one under consideration. From these numbers we calculate the 
effective efficiency of the cut. The differences between data and MC are taken as 
the systematic errors which are added in quadrature. Results are shown in Tables. 
17 11 ^ 17-21 The 4-momentum cut does not appear because its efficiency is close to 
100%. 



Cuts 


MC eff. 


Data eff. 


Systematic Error 


eOverP2 


96.2 


96.6 


0.4 


QED e^e^ 7e'^e^ suppresion 


93.1 


94.2 


1.2 


QED /i+yU^ — > 'jp~^p~ suppresion 


75.3 


77.7 


3.2 


Hadron separation 


96.9 


96.8 


-0.1 


Overall analysis cut systematic error 






±3 



Table 7-1: Efficiencies for T(15') 77r+7r^ for flat signal MC, efficiencies from 

data (p), and the derived systematic error in %. Efficiencies are reported as the 
number of signal events after all cuts divided by the number of signal events with 
all cuts except the one under consideration. Statisticall errors are 0.1% or less. 



For the proton case we don't have a clean sample with high statistics of 
e"^e~ — > events in data. By extension we take the systematic error in this mode 
to be 10%. 

7.1.1 Justification of the DPTHMU Cut 

The reason we prefer to use DPTHMU < 5, instead of the more traditional 
(see j^) DPTHMU < 3 used in CLEO II, is that, for some unknown reason, our 
CLEO III MC has too many vr tracks with 3 < DPTHMU < 5. 



^ A study using or A signals from hadronic environments would have larger 
statistics, but is problematic because the large number of tracks and showers artifi- 
cially decrease the cut efficiency. See j2H] Appendix A. 3 for an example. 
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Cuts 


MC Eff. 


Data Eff. 


Systematic Error 


eOverP2 


98.5 


100 


1.5 


QED e+e^ 'ye~^e~ suppresion 


98.1 


99.3 


1.3 


QED fi'^fi" — > 7/i^/i~ suppresion 


93.5 


96.4 


3.1 


Hadron separation 


88.4 


82.2 


-7.0 


Overall analysis cut systematic error 






±8 



Table 7-2: Efficiencies for T(15') 'jK^K^ for flat signal MC, efficiencies from 

data (0), and the derived systematic error in %. Efficiencies are reported as the 
number of signal events after all cuts divided by the number of signal events with 
all cuts except the one under consideration. Statisticall errors are 0.1% or less. 



To observe this fact we first select a relatively clean sample of pion tracks 
by requiring the event to have a tt+tt^ invariant mass consistent with the p mass, 
and to pass all our analysis cuts except for the cut on DPTHMU on one track. 
For such events we plot the rate as a function of momentum, at which the vr track 
whose DPTHMU cut we released has DPTHMU > 3 and DPTHMU > 5. To 
increase our statistics we do this procedure twice, once for each track, and average 
the fake rate. Figure 17^ shows the results. Clearly the CLEO III MC we are using 
has some problem modeling the DPTHMU < 3 cut. 

To keep the systematic error low we choose a cut at DPTHMU < 5 
{DPTHMU < 3 gives a systematic error of about 20%). This does not change the 
efficiency in data very much, but it increases the efficiency reported by MC, bring- 
ing it closer to reality. The increase in /x fakes after loosening the cut is estimated 
to be low using QED MC (see Figure I7~2|l . 

7.2 Angular Distribution of Signal 

The photon and tracks from the process T(IS') —>■ 7X with X — > h'^h" have 
a different angular distribution than that of flat MC. Examples of possible angular 
distributions are shown in the appendix. 
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File: 7bohr/LJser1/luis/cleo3/ntp/data_50.ntp 



ID 
1 
1 
1 
1 



0.20 



0.10 



3 



0.00 



0.040 - 



0.020 



0.000 



IDB Symb 

4 31 
3 -31 
6 31 

5 -31 



Date/Time 
040525/1408 
040525/1408 
040525/1408 
040525/1408 



Area 
0.2350 
0.4317 
9.6769E-02 
6.051 6E-02 



Mean 
2.758 
2.684 
2.996 
2.792 



1.6 

Track Momenta (GeV) 



3.2 



R.M.S. 
0.8547 
0.8088 
0.7400 
0.6995 




Figure 7-1: Pion faking muon fake rates for a cut on DPTHMU < 3 (top) and 
DPTHMU < 5 (bottom). MC is shown as solid circles while p from data is shown 
as hollow circles. Fake rates are reported relative to all events that pass all cuts, 
except for the DPTHMU cut for one of the tracks. 



Figure 17 31 shows the efficiency in flat MC as a function of cos 6ry and 
cos6h+'^ . Note that the K and p modes are nearly insensitive to the track an- 
gular distribution, while the vr mode is more sensitive to 9^+. This happens because 
of the stronger muon rejection cut in the vr mode. 

We measure the systematic effects of fiat MC efficiency by convoluting each 
plot in Figure 17^ with different possible angular distributions calculated in the 



^ 6^ and Oh+ are the helicity angles of the sequential decay, defined in the ap- 
pendix (see Figure |21). 
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File: */bohr/user1/luis/cleo3/ntp/mcpass2_mumu_gamma_new.ntp 

ID IDB Symb Date/Time Area Mean R.M.S. 

1 3 -31 040518/1402 3.676 2.361 0.9088 

1 30 -32 040518/1402 3.022 2.549 0.8590 




Track Momenta (GeV) 

Figure 7-2: Muon faking pion fake rate (in %) for DPTHMU < 3 (solid squares) 
and DPTHMU < 5 (solid circles). There is a complete overlap in the last bin. 



appendix. Tables 17-3117-41 show the necessary correction factors relative to fiat MC 
efficiency for decays with definite 7 and X helicities and Ax, due to the non-flat 
photon and hadron distributions. We call these factors ^x^Xx^"^^ ^^'^ ^a^a^^^) 
respectively. 

The efficiency of a decay with definite A^, Ax can be obtained using the fiat 
MC efficiency corrected by a factor e'l^Xx ~ ^a^Ax(^) ^ ^a^Ax(^)' general, the 
final state is a mixture of all possible A^, Ax pairs, and the efficiency correction 
factor is, 



s-^^ = cos^ e cos^ $ei^ + sin^ Ge^f + cos^ 9 sin^ ^e(i . (7.1) 
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Figure 7-3: Flat MC efficiency as a function of cos6^ (left column) and as a func- 
tion of cos6h+ (right column). Top row corresponds to the pion mode, middle row 
to the kaon mode and bottom row to the proton mode. 



The fits in Figures f4-8l and 14-lUI and I5-6H5-8I measure the pair (O, $). These 
values are summarized in Table 17-51 where they are used to obtain e"^^ for each 
mode. 

The pair (B, $) carry an error which is a source of systematic uncertainty. 
We calculate this systematic uncertainty by inspecting the differences in efficiency 
when (0, $) move away from the value which gives the minimum chi-squared, Xmin^ 
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Mode 


aio or ai2 


ail 




= 


e = - 


TT 


0.910 


1.18 


K 


0.925 


1.15 


P 


0.922 


1.16 



Table 7-3: e^^^x^^)' efficienciy correction with respect to flat MC factors due to 
the non-flat angular distribution of the photon in T(IS') 'jX for different X spin 
values (Jx) and 7, X hehcities (A^, Ax)- 



Mode 


Jx 


aio 


an 


ai2 






= 0, $ = 


6 = - 

^ 2 


e = 0, $ = f 







1.00 








1 


0.785 


1.11 




X TT + TT- 


2 


0.829 


0.934 


1.15 




3 


0.898 


0.853 


1.03 




4 


0.899 


0.863 


0.922 







1.00 








1 


0.944 


1.03 




X K+K- 


2 


0.950 


0.988 


1.04 




3 


0.957 


0.971 


1.01 




4 


0.965 


0.962 


0.992 



Table 7-4: e^^^x^^)' sfficienciy correction factors with respect to non-flat MC due 
to the track angular distribution in T(IS') 7X, with X — » h'^h^ for different X 
spin values (Jx) and 7, X heli cities (A^, Ax)- 

under the condition < Xmin + 1 ^ - For Jx > 1 both (B, $) are free to move, 
defining a surface in the (0, $) plane. These surfaces are shown in Figure 17-^ for 
the /2(1270), /4(2050), and /si 1525). 

At this point we can check whether (B, $) depend on the mass of the de- 
cay. We split the /2(1270) and the /2(1525) into a high mass and a low mass 
region. The plots of the error surfaces of the measured (0, $) show no significant 
separation for the different mass regions (see Figures 17-51 and l7-6j) . 



^ In a two dimensional linear problem such a set of points defines the surface of 
the standard error ellipse. 
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Results for the correction factor and its systematic error are shown in Ta- 
ble Ol 

Upper hmits on T{1S) ^ 7/2(2220), /2(2220) h+h- are changed to 
include the angular distribution's effect on efficiency. Since we can't measure 
the helicity amplitudes in this case, we choose the worst possible case where the 
corrected efficiency is lowest. This always corresponds to = $ = 0. Results are in 
Table 01 
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Figure 7-4: Surfaces in the (0, $) plane used to determine the systematic uncer- 
tainty in the efficiency correction factor for the modes with /2(1270), /4(2050), 
and /2(1525). The /4(2050) surface may seem large, but when drawn in spherical 
coordinates it is a small "north pole cap" . 
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Figure 7-5: Measured (G, $) surfaces for the /2(1270) high and low mass regions. 
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Figure 7-6: Measured (G, $) surfaces for the /2(1525) high and low mass regions. 



Table 7-5: Measured (O, $), calculated e^^ , efficiency correction factor interval when (O, $) move away from their measured 
value under the condition that stay within one unit of its minimum value, and the systematic uncertainty on the correc- 
tion factor. 



Mode 


Jx 


e 


$ 




Interval 


Sys. Err. (%) 


7/o(980) 











0.91 


0.91 - 


0.91 





7/o(980) 


1 







0.71 


0.71 - 


0.80 


+13 


7/2(1270) 


2 


'-'•'-"-'-0.00 


'"'•'J '"'-0.10 


0.78 


0.76 - 


0.80 


±3 


7/4(2050) 


4 


1 cry+O.OO 
' -0.34 


0.2 ±0.7 


1.02 


1.00 - 


1.02 


-2 


7/2(1525) 


2 


U-UU-0.00 


'"'•'J'J-0.16 


0.90 


0.89 - 


0.91 


±1 


7/o(irio) 

7X+K-, 2 GeV/c^ < m(K+K-) < 3 GeV/c^ 
7PP, 2 GeV/c^ < m(pp) < 3 GeV/c^ 











0.93 


0.93- 


0.93 





1 


'"'•-'-'-'-0.16 





0.88 


0.87- 


0.91 


+3 
-1 


1 


U.UU_Q QQ 





0.92 


0.92 - 


0.95 


+3 



Table 7-6: Worst case efficiency correction factors and their effect on upper limits of T(15') — > 7/j(2220), /j(2220) — > h^h 
decays. 



Mode 


Worst case correction 


Upper limit increase (%) 


/j(2220) ^ TT+TT- 


0.75 


+33 


/j(2220) ^ K+K- 


0.88 


+14 


/j(2220)^pp 


0.92 


+9 
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7.3 Different Hadronic Fake Rates Between IS and 4S 

Differences in tlie fake rates of liadronic events between tlie IS" and 45" 
data are anotlier potential source of systematic errors wlien doing a continuum 
subtraction. 

In principle this is a second order effect since the contamination is proportional 
to the difference of two small numbers. However since we have such a large number 
of p and events in the underlying IS continuum, we should quantify any possible 
signal contamination. Section 3.2.2 outlines the three cases we need to worry 
about; p events contaminating the K^K'^ invariant mass plot, p contaminating the 
pp invariant mass plot, and events contaminating the pp invariant mass plot. 

It is reasonable to assume that the IS and 4S fake rates stay within 50% 
of each other (see Section 3.3.2 for some examples from MC measurements). 
With this assumption, taking np ^ 20000 and ^ 2000 in the IS data, and 
the fake rates measured in data (see Table l3-(i|l we can calculate the systematic 
uncertainties. We expect ±27 events from p to contaminate the K^K~ invariant 
mass plot. However, these events can be ignored as source of systemtic uncertnty 
because 90% of them fall in the 1 — 1.5 GeV/c^ region, below the /2(1525) peak. 
We also expect ±7 from p and ±3 events from contaminating the pp invariant 
mass plot. This represents a ±8% systematic error for the T(IS') — > 7pp mode. 

7.4 Other Systematic Sources 

Besides the systematic uncertainties described above, we also add a 2% 
systematic effect from track finding (1% per track) and 5% from the number of 
T{IS). 

We find no evidence for non-resonant hadron pairs that could interfere with 
the resonances. There is a possible interference between /2(1275) and /2(1525) 
because they both are in the J = 2 state. From the /2(1275) measurement in the 
TT+TT^ mode we expect 3 events from this source in the /2(1525) mass region 1.45 — 



103 

1.6 GeV/c^. In the most extreme cases we have, oc 1 1^/^(1525) I P + 1 1^/2(1270) I P =t 
2| 1^/2(1525) 1 1 1 1^/2(1270) 1 1) where ^^^(1525) and A/2(i27o) are the amphtudes associated 
with /2(1525) and /2(1270) respectively. From the data we have ~ 300, and 
using the fact that 1 1^/2(1270) I P contributes with 3 events, the contribution from 
11^/' (1525) 11^ can range from 400 to 200 events. This represents a 30% variation from 
the central value. Since this is the extreme interference case, we take half of this 
value, 15%, as a reasonable systematic uncertainty from possible interference in the 
/2(1525). 

The possible decay T(IS') — »• pTT could contaminate our vr"*" tt" invariant mass 
plot, since a vr'' with a momentum greater than 4 GeV/c^ looks very much like a 
photon. Using phase space MC and isospin symmetry, we estimate that about 40% 
of possible T(IS') pvr have such a 7r°. The latest measured upper limit j211 
for this decay mode is 1.9 x 10~^. This gives an upper limit of 16 events in the 
efficiency corrected vr^ vr^ invariant mass plot, of which 10 events fall in the p mass 
region, behind the /o(980) mass region. The remaining potential 6 events could 
contaminate the rest of the invariant mass region. None of the signal regions is 
affected by more than 1%. 

To test the robustness of the statistical fits, we redo them with all masses and 
widths of possible resonances required to stay within la of their nominal PDG 
values [^. The results from fitting the data in this way are within the statistical 
errors when compared to the fit results summarized in Table 14-11 and 15 II We do 
not observe a systematic difference between these two fitting techniques. 

The branching ratios for /2(1270) tt+tt-, /4(2050) tt+tt^, and /2(1525) ^ 
K^K^ are taken from the PDG ^E] and contribute to the systematic uncertainty 
with 1^3, ±9, and ±3 percent respectively. 
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7.5 Overall Systematic Uncertainties 

We combine the systematic errors of each section in quadrature. Tables 17-71 - 
17-81 summarize the systematic errors. 

Table 7-7: Systematic errors expressed as % from the different sources described 
in the text. Not shown in the table, but included in the total, are a systematic 
uncertainty of 5% from the number of T(IS'), a 2% systematic uncertanty from 
MC tracking in all the modes, a 15% systematic uncertainty in the 7/2(1525) from 
possible interference with 7/2(1270), an 8% systematic uncertainty in the 'ypp 
mode from possible hadronic contamination, and the systematic uncertanties in the 
/2(1270), /4(2050), and /2(1525) hadronic branching fractions. 



Mode 


Analysis 
Cuts 


Angular 
Distribution 


Total 


7/o(980) 
7/2(1270) 
7/4(2050) 


±3 
±3 
±3 




±3 
±2 


±6 
±7 
±7 


7/2(1525) 
7/o(1710) 
-fK+K-{2-3 GeV/c2) 


±8 
±8 
±8 


±1 



+1 

-3 


±20 
±10 
±10 


7pp(2 - 3 GeV/c2) 


±10 


-3 


±14 



Table 7-8: Increase in the upper limits of T(15) 7/j(2220) for different /j(2220) 
decay modes due to the possible systematic effects described in the text added in 
quadrature. Not shown in the table, but included in the total, are a 5% contribu- 
tion from the number of T(IS') and of 2% contribution from MC tracking. 



Mode 


Analysis 


Angular 


Total 




Cuts 


Distribution 




/j(2220) TT + TT- 


±3% 


±33% 


±34% 


/j(2220) ^ K+K- 


±8% 


±14% 


±17% 


fj{222Q)^pP 


±10% 


±9% 


±14% 



CHAPTER 8 
RESULTS AND CONCLUSION 

We report on a new search for two-body radiative T(IS') decays. We place 
stringent limits on the production of the /j(2220) particle in the pion and kaon 
modes, and a less stringent limit in the proton mode where some excess of events is 
observed in the region of interest. 

In the decay channel T(IS') 77r+7r^ we find clear evidence for the resonance 
/2(1270) and measure a branching fraction of (13.3 ± 1.4) x 10^^, which is 
consistent with the earlier CLEO measurement of (TA'tfl) x 10~^ ^H]- The angular 
distributions of the photon and tracks strongly indicate that the hadron pairs we 
assign to the /2(1270) are indeed in a J = 2 state, and that nature prefers to 
produce the /2(1270) with hehcity. In contrast, for J/ip 7/2(1270) it was 
found that the /2(1270) is produced at equal rates with both helicity and 
helicity 1, and at a rate consistent with for helicity 2, but more recently ^20] 
measured helicity dominance for this same mode. 

There is a barely significant (4.3o") excess of events in the /o(980) —>■ tt+tt^ 
invariant mass region. The angular distributions of the photon and tracks indicate 
that this excess of events is in a J = 1 state, rather than J = which would be the 
case if the excess were due to /o(980) decays. It may be that this excess is due to 
non-resonant 77r+7r^ production. We conclude that more data is needed to resolve 
this situation. 

We find weak evidence for the production of the resonance /4(2050) in the 
T(IS') — > 77r+7r^ decay channel. 

In the decay channel T(IS') '~fK^K~ we find strong evidence for the 
production of the resonance /2(1525) and weak evidence for the production of 
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/o(1710). The photon and track angular distributions show that the two tracks 
attributed to the /2(1525) are indeed in a J = 2 state and that the /2(1525) 
is produced mostly with helicity 0. This is consistent with Kramer's theoretical 
prediction and similar measurents done in the J/ip system |32j . 

There is also evidence of an excess of events in the 2 — 3 GeV/c^ K^K^ 
invariant mass region, which we cannot attribute to any known resonances. 

Finally, we also find evidence of an excess of events in the channel T(IS') — ^ 

which we cannot attribute to any known resonances. We find no evidence 
of an enhancement near pp threshold. Such an enhancement was recently re- 
ported by [Sn] in the J/ip system, and is currently being interpreted (See 
e.g. [SI, M)- 

When comparing the radiative T(IS') decays studied in this analysis to ra- 
diative J/ip decays we observe a suppression ratio of 0.09 that is in reasonable 
agreement with with the naive expectation of 0.025 obtained from scaling argu- 
ments. 

Table 18-11 summarizes our results. 

While we find hadronic resonances at a level consistent with the estimate 
of Section 1.2, we do not observe any glueball states. This is not the result we 
would expect from a naive interpretation of QCD, as the two gluons from a T(1S) 
radiative decay can form a bound state directly and need to go through at least 
two strong interaction vertexes to form a hadron. 

Internal consistency in QCD predicts the existence of glueballs. Even though 
theoretical calculations of the glueball mass spectrum exist, a clear glueball 
observation has not been made yet. We conclude that more glueball searches are 
necessary to clarify this situation. 



Table 8-1: Final measured branching ratios, measured branching ratios reported relative to J/ if) branching ratios, and sta- 
tistical significance for each decay channel. The measured branching ratios have been corrected by the factors calculated in 

Appendix A2. For the branching ratios the first uncertainty is statistical and the second is systematic. 



Channel 


Branching Fraction x(10 5) 


T(15) 


Significance 


7/o(980), /o(980) ^ TT+TT- 
7/2(1270) 
7/4(2050) 


i.8i^:«±o.i 

13.3 ±1.0 ±0.9 
3.5 ± 1.3 ±0.2 


0.09 ±0.01 
0.019 ±0.006 


4.3(7 

> U(x 
2.6(7 


7/2(1525) 
7/o(1710), /o(1710) ^ K+K- 
-fK+K- 


4.3+i:^±0.9 
0.38 ±0.16 ±0.04 
1.14 ±0.14 ±0.11 


007+° °°^ 


> 14(7 

3.2(7 

9.1(7 


rypp 


0.41 ±0.08 ±0.10 


0.011 ±0.005 


4.8(7 


7/j(2220),/j(2220)^7r+7r- 
7/j(2220), /7(2220)^X+ir- 
7/J(2220),/J(2220)^pp 


< 0.08 

< 0.06 

< 0.11 







APPENDIX 

HELICITY FORMALISM FOR TWO BODY DECAYS 
Helicity formalism can be used to obtain the angular distribution of a decay 
process. In this appendix we restrict ourselves to the case of two body decays. 
The situation is as follows: particle A with spin J^, and z-axis spin projection 
decays to two particles, A ^ B + C with definite helicity and Ac- Let the 
direction of the of the decay be n{6,(j)) = n{9B,(pB) = Pb = —Pc (see Figure H]). 
The amplitude for this decay is, 

A^'X = < H\b\c\H\JaMa > . (1) 

We have little chance of directly calculating ^ But we can exploit the fact that 
H conserves angular momentum by considering the basis {< jmXBXc\}, where j 
is the total angular momentum of B + C and m its projection along the z axis^ . 
Inserting this basis into Equation ^ 

jm 

= X < 9(l)\B\c\jmXB\c > 5jj^5mMA^XB>^c 
jm 

= < e(t)\B\c\JAMA\B\c > AxgXc- 



^ The rotational invariance of helicities is a necessary condition to construct 
{< jmAeAcI} as a basis for the two-particle state 5 + C. Rotational invariance 
is one of two key ideas in helicity formalism. The other key idea, which is boost 
invariance, is useful in sequential decays. 
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Figure 1: B + C. 

The above equation contains two factorized terms. The first term is indepen- 
dent of the interaction H and contains the angular distribution information we are 
interested in. The second term, Ax^x^, is usually called the "helicity coupling am- 
plitude" . It is independent of the angular distribution, and contains the physics of 
the decay (in particular it depends on and A^, but not on Ma)- In the helicity 
formalism, the helicity coupling amplitudes are unknown parameters. 

For strong and electromagnetic decays there is another invariance of H we 
can exploit; parity. For these kind of interactions it can be shown by inserting the 
parity operator that, 

A_x,-xc = {-iy^-'''-'''VAVBVcAx,Xa (2) 

where tjatIbVc is the product of the parity of particles A, B and C. This useful 
property of AxgXc reduces the number of unknown parameters in the angular 
distribution of strong and electromagnetic decays. 

The next step is to derive explicit formulae for < 94>XbXc\JaMaXbXc > for 
the B + C system. Such derivation, which is out of the scope of this appendix, will 
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be omitted here^ . The result can be expressed in terms of the Wigner functions, 



< e<pxBXc\JAMAXBXc >= (s) 

where A = — Ac and, 

Di\,M 0, -0) = e-^^^^--^M'-^^,{e). (4) 

The functions c?]^^ ^(^) have real values, and some of them are given in the particle 
data book. 

We have finally arrived to the following important relation. 



^"^c = \r-^i>Mi.A,^^c- (6) 



All helicity distribution calculations are reduced to handling these basic ampli- 
tudes. Here are some of their very important properties, 

A^X = ^ \^b-Xc\>Ja, (6) 



A-^{^ 11 = 11 Af-\ II, (7) 

A^Ac II II A^Ac II' V y 



27r / / / 



(^^^^XnXc\'^x'^ - \x--X') II ^>^bXc II • (9) 

Property El is a reflection of the physical fact that the projection of the total 
angular momentum along a particular axis cannot be greater than the total angular 
momentum itself. 



^ The interested reader can find this derivation in j^Tj 
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II '^'xbXc 11^ angular distribution of a final state where and the 

helicities of B and C are known (see Equation Q). However, it is usually the case 
that MyiAfiAc are not known. In such cases we can only measure the overall 
angular distribution, which is simply the average over possible initial states of the 
incoherent sum of final states, 

E E ^«.JI<A.II^ (10) 

where Pma is the weight of the initial state Ma and dQs = d{cos{9B))d(f)B- 

Next, consider a sequential decay, A ^ B + C with C D + E (see Figure 
El)- The rotational and boost invariance of helicities makes it easy to extend the 
helicity formalism to this case. To do so we consider a new reference frame, O', 
at rest with C, such that z' = n{6B,(f)B)^ ■ All particle helicities remain the same 
in O', and in particular Mz{C) = —Xc (the negative sign arises from the fact 
that n{6B,(f)B) is the direction in which B is traveling). It is now straightforward 
to calculate the intermediate amplitude C^^J^ at O'. Of course ^aI'ac ^^^^^ 
calculated at A's original rest frame. Since Ac is an intermediate state, the helicity 
amplitude for the final state A^AdA^; is, 

AC^\ . =yA^'\ Cl'^'l (11) 

Ac 

where we have coherently summed over the intermediate state C . If M^A^AdA^ 
are not known, 

^ P r\ P^^^ 

p "-^AflApAfi _ p II ,pMA ||2 /-iQN 



d^Bd^D d^Bd^D 

Ma^b^d^e Ma^b^d^e 



^ This condition does not unlikely determine O'. We can be more precise by 
defining O' = L{'pc)R{4>B, ^b, 0)0, where L is the Lorenz Boost operator, R is the 
rotation operator in terms of the Euler angles, and O is A's original rest frame. 
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Figure 2: A ^ B + C and C —>■ D + E. The drawing is peculiar in that C's daugh- 
ters are drawn as seen in C's rest frame (back-to-back). 

The sequential decay amplitudes have the following important properties, 



AC 



-Ma 



AC 



Ma II 



(13) 



[\<pj ^d<p, II AClX,^ r= An^J2 II A^X 111 C^S.. f • (14) 
^0 Jo \ 

Ac 

This kind of logic can be extended to a general two-body decay chain. 
Example: T{1S) 7/2(1270)(/2(1270) ^ tt+tt-) 

From our previous discussion we can immediately write (see Equation ITTj) . 

^iK+K- Z-^ -^7^/2(1270) ^^+K- 

A/2(1270) 

In this case, A^- = = 0, = ±1. In CLEO Afx(i5) = =tli there is no 
^^T(i5) = since at the T(IS') energy electrons couple only to positrons of the 
opposite helicity, this way Pi = P_i = 1 and Pq = is an excellent approximation.. 
There are therefore four amplitudes we must consider, w4CJoo5 -^^roo' -^^-loo 

■^^-100- 
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The overall angular distribution is (see Equation IT^ . 

— 1 II /Ifl ||2 ,1 II Arl ||2 ,1 II A^-l ||2 ,1 II Ar-^ ||2 

— II AC^ IP -I- II AC^ IP 

— II ^tL^ioo II T II >lL-_iQQ II 

where we have used property [121 The two relevant amplitudes are, 



•^'-'100 ~ ■^lO'-'OO ^ ■^ll'-'OO ^ •^12'-'00 

Afi — /ll I /ll -I- /l^ C 

•^'-'-100 ~ ■^-lO'-'OO '^-l-l'-'OO ■^-1-2'-' 



(15) 



2 

00 



where we have used the fact that — A/2(i27o)| < Jr{is) = 1 to sum over all 
possible A/2(i27o)- 

All that remains is to substitute, 

II Afl ||2_ II /ll (-0 ||2 , II Al n-l ||2 , || aI n-2 ||2 , 
II -^'-'100 II ~ II ■^lO'-'OO II ^ II '^ll'-'OO II ^ II ■^12'-'00 II ^ 

i r° A*^c*^^ -)- A^ r° /i*if *-2 I /ll f-i /i*i(^*-2 I „ „ 1 

I'^lO'-'OO'^ll'-'OO •^10'-'00'^12'-'00 ■^ll'-'OO '^12'-'00 L.C.J- 

II /1(-1 ||2_ II Al ||2 , II Al ^1 ||2 , II Al ^2 ||2 , 

II '^'-'-100 II ~ II '^-lO'-'OO II ^ II '^-l-l'-'OO II ^ II '^-l-2'-'00 II ^ 

I A^ r° A*^ C*^ -\- A^ r° A*^ C*"^ A- A^ A*^ C*'^^rr\ 

I'^-lO'-'OO'^-l-l'-'OO ^ '^-10'-'00'^-l-2'-'00 ^ '^-l-l'-'00'^-l-2'-'00 ^ c-L-J 

(16) 

into Equation UHl 

Doing so will give us the full angular distribution of the decay chain. However, 
using property |Hl all the cross terms inside {...} vanish if we integrate over 0^+ or 
(f)^. For the time being we integrate over both 0^+ and 0^. We will return to the 
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terms inside {...} later. 

r-2-n r-2-K 

Jo ^'^ Jo '^'^ dil^dilT^+ 

^Aiw IP +11 AUo in iidp + 

(II Al, r + II A\_, IP) II Cl, IP + 
(II Al, IP + II A\_, IP) II Cl, IP] 



= f [II ipii Coo IP {d\Ao,r + d\_,{9,r)di{9^,r + 

II All IPII Cqo IP ((i}o(^7)^ + '^lo(^7)^)'^lo(^7r+)^ + 

II Ai2 IPII Coo IP {dU{0,Y + dl_i{9,y)dl,{9^+Y] 



= ^ II Coo IP [II Ao IP ^(l + cos2^,)(^cos2^.+ - + 

II All IP sin^ ^-yiy 2 ^i^^7r+ COS^7r+)^ + 

II Ai2 |pi(l + cos2^-,)(^sin^^,+ )]. 
Z o 

At this point it is convenient to define the normalized helicity amplitudes, 

^10 

Oio = 



On = 



ai2 = 



VII ^10 


IP + II Ai IP + 1 


1 Ai2 IP 




An 




VII ^10 


IP + II IP + 1 


1 Al2 IP 




Al2 





Uo IP + II All IP + II Ai2 IP 



which satisfy. 



oio IP + II On IP + II 0,12 IP= 1- 
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The normalized probability distribution, ^^(cose )d{cose +) ' obtained by 

dividing the probability distribution by the branching fraction of the sequential 
decay, 



P 



2 II Coo f (II Ao f + II An f + II ^2 



Coo iHl ^10 |P + II Coo iHl All iP 



Coo iHl Ai2 Ip 



+ 



Coo iHl ^-10 IP 



Cqo IPII v4_i_i IP 



Cqg IPII 11^ 



This way. 



d{cos ^^)(i(cos6'^4 



flio 



;i + cos^ ^^)-(3 cos^ 9^+ - If + 



Ii2 ^ '2/1 / • \ 2 

ail II ^^1^ c^^ — (sint^7,+ cosc/^+) + 

ii2 ^ /-I 2 n \ • 4 /I 

ai2 -(1 + cos 9^)—sm 9^+. 
8 lo 



It is useful to project the angular distribution on 9j and 6'^ 



dU 



d cos 9^ 



(^10 IP -(3cos^6'^+ — 1)^+ II ail IP -r(si^^7r+ cos^^^+)^+ 
8 4 

II 2 -'-5 • 4/1 

ai2 77 sm 6/^+ 
lo 



(18) 



dU 



d cos 6*^ 



3 3 
flio f + II ai2 f )-(l + cos2^^)+ II ail f -sin^^^. (19) 



Let us now return to the terms inside {...} in EquatiorlKil It turns out that for 
these terms a substantial simphfication occurs if we integrate over cos6'7r+. After 
some algebra, we arrive at the following expression. 



= ^(1 - -^{aioa*2 + atoai2} cos|2(0^ - 0^+)]). (20) 
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Example: T{1S) 7/4(2050) (/4(2050) ^ tt+tt-) 

For the case /4(2050) the calculation is very similar to the /2(1270) case 
discussed above. The allowed helicities for /4(2050) are the same as the allowed 
helicities for /2(1270). This makes the angular distribution calculation very similar 
to the /2(1270) case, the only difference being that the functions 2(i270)0 
replaced hy di r,. 

The photon's angular distribution comes out to be the same. The result for 
the pion's normalized angular distribution is, 

'^^ " aio f —(35 cos^ - 30 cos^ + 3)^ + 



(icos6'7r+ 128 



flu -^(sin6'^+ cos6'^+)^(7cos^ 6'^+ —3)^ + 
ai2 IP ^sin'^6'^+(7cos^6'^+ - 1)^. 



(21) 

Fit of the Helicity Angular Distribution Obtained from Data 

In practice, a particular bin in the helicity angular distribution obtained from 
data has bins with dN = NdU events, where is the total number of events 
(which is allowed to float during the fit), and dU is the portion of the normalized 
probability assigned to the bin by the helicity formalism prediction. 

To incorporate Equation El directly into the fit we make the following change 
of variables, 

II oio 11= cos 6 cos $; || ai2 ||= cos 6 sin $; || an ||= sin 9. (22) 

Because of this definition, both O and $ should be required to be within the 
interval (0, ^) during the fit. The new variables have the convenient property that 
for J=0, $ = = 0, and for J = 1, $=0. In general there are 1, 2 and 3 degrees 
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of freedom for the fits with J = 0, J = 1, and J > 1 respectively. One degree 
of freedom always corresponds to the number of events, A^, and the other two 
correspond to G and $, which determine the magnitude of the normalized helicity 
amplitudes through Equation |221 • 

When fitting, it is preferable to choose the full 2D helicity angular distribution 
which is a function of 9.^ and 9t^+ instead of its one dimensional projections. The 
2D choice is better not only because the 2D distribution carries the maximum 
amount of information, but because in general the efficiency depends on both the 
pion and photon helicity angles, which can lead to complications when efficiency 
correcting the one-dimensional projections. Building on the two previous examples 
the 2D helicity angular distributions for J = 0, 1, 2, 3, 4 are (for compactness we 
define x = cos and y — cos d^)^ 



^^iV,.,(cos^e-(l + ,^)-x^ + 

sin^e^(l-y^)^(l-x^)) 



(24) 



Obtaining the error of the normalized helicity amplitudes requires some matrix 
multiplication. For compactness let us define a = | ai2 | and cu = ( ^ ) • The error 




matrix for the normalized helicity amplitudes is, 

< 5a5a >= — < 600600 > — 
ouj ouj 

where < 6oj6uj^ > is the error matrix obtained from the fit. 
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^ = Nj=,i cos' e cos^ $ ^(1 + y'f-iSx' - 1)'+ 
axay 8 8 

sin^e ^(i_^2)15^2(^_^2)^ (25) 

4 '4 

cos^esin^ $^(1 + 2/^)^(1-0:2)2) 
8 lb 



f^^iv..3(cos2ecos2.^(i + ,2)r^5,3_3,).^ 

sin^e ^(1- ^2)21(1 _^2)(5^2_ 1)2+ (26) 

cos^ e sin^ $ -(1 + y')—x\l - x'f) 



^ = iV,=4( cos2 e cos2 $ ^(1 + 2/^)^(350;^ - 30x' + 3)^+ 
axay 8 128 

sin^ e ^(1 - y')^x\l - x'){7x' - 3)2+ (27) 

COs2esin2$ ^(1 + ^2)45 2)2(7^2 _ i)2)_ 

8 64 

In this analysis, because of our low statistics (and by low statistics we mean 
that there are numerous bins with < 1 events in the 2D plots), we project Equa- 
tions EBIEZI on X = cos^TT, and y = cos^^. We also fold each distribution around its 
symmetry axis x = and y = to increase each bin's statistics. Finally, we do a 
simultaneous fit to both distributions and obtain N, 9, and $. 

Admittedly, this procedure has the problems described previously when 2D fits 
where advocated, which can give rise to systematic effects in the measured helicity 
amplitudes. In essence, we need to use a unitary weight function proportional 
to the efficiency distribution (see Figure l7-3p in the integrals used to make the 
projections. This affects the relative strength of each helicity amplitude and the 
overall efficiency, but does not change each projected shape. However, the helicity 
amplitudes themselves are only used to calculate systematic effects on the efficiency 



119 

of resonances, so by simultaneously fitting the unweighted ID projections we are 
just ignoring the systematics on the systematics. 

Finally, let us consider the 0^ — 0,^+ distribution for J = 2. Equation EOl 
becomes, 

."^^"t , = %^(1 - ^{cos2esin2<|.cos A}cos[2(0, - M]) (28) 

"(07-071+) 277 y/Q 

where the new degree of freedom A is the relative phase between aio and 020- The 
deviation from the flat distribution is characterized by 

= ^{cos2esin2$cosA} < (29) 

The bigger this factor is, the easier it is to observe in data. Figure El shows the 
continuum subtracted 0^ — 0,^+ distribution for events near the /2(1270) resonance. 
The factor R obtained from fitting the distribution is non-significant. It is a bit 
tedious but straightforward to calculate the 0^ — 0,^+ distribution for other J values. 
We will not present those results here. They are not very interesting because we 
need large statistics to observe deviations from the fiat distribution, and there is no 
effect on the efficiency from non-flat 0^ — 0,^+ distributions. 
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MINUITx'Fitto Plot 9&3 

UlS^Gamma h h. pipiphi axis 

File: Generated internally 

Plot Area Total/Fit 779.39 / 779.39 

Func Area Total/Fit 769.52 / 769.52 

5.2 for 10- 2d.o.f., 

Errors Parabolic 
Function 1: COMIS Function XMNCMJ 
N 769.49 ± 56.78 

R -6.87430E-02 ± 0.1062 

150 
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Fit Status 3 
E.D.M. 1.509E-06 
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Figure 3: Continuum subtracted 0^ — 07r+ distribution for events near the /2(1270) 
resonance with a fit to Equation |2H1 overlayed. 
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Example: Continuum Events. 

Our previous examples tackled radiative decays of the T(IS'). Now let us 
consider the continuum "radiative beam" process e+e^ —>■ 77r+7r^. The helicity 
formalism can be applied at the CM of the tt+tt^ vertex, where the initial state has 

= 1~ . We have 1^ 0^0^ at the tt+tt^ vertex and there is only one helicity 
amplitude we need to consider, 

^00 = \/J^to^o (30) 
which after a couple simple steps leads to, 

iV^sin2^.+ . (31) 



d cos 9t^+ 4 

This result is also valid for the kaon case. 

On the other hand, the photon's angular distribution cannot be calculated 
using two-body helicity formalism. A good approximation is (see |i2Bj Appendix B), 



dP 1 + cos2 



d cos 6 J 1 — cos^ 6y 
which is not valid at low 6 angles (6' — 0, ^ vr) because of the effects on the 
approximation of the neglected electron mass. 

Figure ID shows these angular distributions in 4S "radiative beam" events. 
The agreement between theory and experiment is excellent. This is reassuring, 
because in the main analysis we rely on our angular distribution measurements to 
get the correct efficiency in data and to identify the J value and helicity amplitude 
distribution of each resonance. 
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cos(9) cos(9) 



Figure 4: Helicity angle distribution for "radiative beam" continuum events 
e'^e^ — > 7p and e^e~ — ^ 70 in the 4S data. Top left (right) plot shows in gray 
the events selected as p tt+tt^ (0 K^K~). Middle left (right) plot shows 

the 7r+7r~ {K^K~) helicity angle distribution for the selected events, and has the 
helicity formalism prediction, Equation |^ overlaid. Bottom left (right) plot shows 
the 7 angular distribution for the photon accompanying the p (0), and has the first 
order QED prediction. Equation 1221 overlayed. 
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